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The  work  reported  in  this  document  was  undertaken  as  part  of  the  examination  of  the 
soundness  of  the  current  standards  for  the  spacing  between  parallel  aircraft  routes 
and  the  enhancement  of  analytical  methods  to  evaluate  future  standards.  This 
interim  report  describes  work  completed  to  date  on  the  Conflict  Monitoring  Parallel 
Route  Spacing  Analysis.  This  analysis  assesses  the  potential  for  collision  and  the 
controller  workload  associated  with  aircraft  flying  on  same  direction  parallel 
routes.  To  assess  the  potential  for  collision  the  analysis  considers  a  conflict 
alert  function  similar  to  that  employed  in  the  National  Airspace  System.  The 
conflict  alert  function  detects  pairs  of  aircraft  which  are  projected  to  violate 
the  radar  separation  standard  within  a  given  time  period.  In  the  analysis  the 
event  of  a  conflict  alert  is  followed  by  a  probabilistic  delay  and  a  resolution 
maneuver  characterized  by  a  randomly  chosen  horizontal  turn  rate.  The  controller 
intervention  rate  is  estimated  by  using  a  simulation.  Actual  aircraft  tracks  were 
sampled  from  the  FAA  data  base  which  supports  this  activity.  These  tracks  are 
initiated  on  the  routes  based  on  randomly  chosen  sector  entry  times  which  reflect 
Che  level  of  route  loading.  For  both  the  potential  for  collision  and  the 
intervention  rate,  trial  results  based  on  a  subset  of  Che  FAA  data  are  given. 
Further  analysis  is  required  to  investigate  opposite  direction  and  transitioning 
traffic.  In  conjunction  with  this  work,  che  reliability  of  the  surveillance  ar.d 
control  systems  has  to  be  addressed  as  well  as  other  performance  measures. 
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INTRODUCTION 


In  November  1976,  the  FAA  Associate  Administrator  for  Air 

Traffic  and  Airway  Facilities  requested  assistance  from  the 
Associate  Administrator  for  Engineering  and  Development  in 
certain  analytical  activities  relating  to  air  traffic 
separation. ( i )  In  part,  that  request  asked  for  an  examination 
of  the  soundness  of  the  current  standards  for  the  horizontal 
separation  of  aircraft  in  the  continental  U.S.  The  request 
also  called  for  an  enhancement  of  analytical  methods  for  the 
operational  evaluation  of  future  standards. 

The  response  to  that  request  is  a  program  within  the  FAA's 
Office  of  Systems  Engineering  Management  (AEM-100)  to  study 
VOR-defined  air  route  separation.  This  study's  initial  goal  is 
to  develop  an  understanding  of  the  relationship  of  safe  route 
spacing  to  system  performance  on  the  high  altitude  CONUS  en 
route  airways.  The  system  consists  of  both  the  airborne  and 
ground  elements  of  navigation  and  air  traffic  control.  After 
the  safety/per formance  relationship  is  better  understood, 
improved  specifications  of  navigation  and  control  system 

performance  needed  to  support  specific  route  spacings  can  be 

developed. 

The  FAA's  VOR-defined  air  route  separation  program  is  based  on  a 
data  collection  followed  by  modelling  and  analytical 
activities.  The  precursor  to  the  data  collection  was  a  mini 
data  collection  in  1975  done  by  MITRE  with  support  from  ANA-220 
at  NAFEC. (2)  From  this  experience,  MITRE  wrote  the 
specifications  for  the  main  data  collection. (2, 3)  The  main 

data  collection  was  planned  and  conducted  by  NAFEC  (ANA-220) 
from  September  1977  to  April  1978.(4)  At  the  present  time, 
NAFEC  is  reducing  the  data  and  compiling  the  data  base. 

Concurrent  with  the  data  collection,  there  are  several  analyses 
being  performed  which  address  the  relationship  of  navigation  and 
ATC  system  performance  to  safety  of  operations  on  the  VOR  route 
system.  These  analyses  address  the  potential  for  collision 
between  aircraft  assigned  to  different  routes  under  various 
conditions.  NAFEC's  analysis  addresses  the  potential  for 
collision  between  aircraft  assigned  to  parallel  routes  under  the 
assumption  that  there  is  no  radar  being  used  to  separate  the 
aircraft. (5)  There  is  also  an  effort  at  Princeton  University 
to  address  the  potential  for  collision  of  aircraft  on 
intersecting  routes  where  no  radar  coverage  is  available. (6) 
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The  first  volume  of  this  report  describes  MITRE's  Conflict 
Monitoring  Analysis,  This  analysis  addresses  the  potential  for 
collision  and  the  controller  intervention  rate  for  aircraft 
assigned  to  same  direction  parallel  routes  when  the  controller 
monitors  aircraft  movements  with  radar  surveillance.  The 
appendices  in  this  volume  present  the  details  of  the  analysis 
performed  to  estimate  the  probability  of  horizontal  overlap  and 
controller  intervention  rate. 

Appendix  A  describes  the  conflict  region  boundary.  This 
boundary  demarcates  those  pairs  of  aircraft  which  are  projected 
to  be  in  conflict  from  those  pairs  of  aircraft  which  are  not 
projected  to  be  in  conflict.  Appendix  B  describes  the 
estimation  of  the  probability  of  horizontal  overlap  given  that 
the  aircraft  pair  is  on  the  conflict  region  boundary.  In  order 
to  make  this  estimate,  one  needs  to  be  able  to  find  the 
probability  of  being  observed  within  the  conflict  region  when 
there  are  uncertainties  in  the  position  and  velocity  estimates 
from  the  tracker.  The  estimation  of  this  probability  is  given 
in  Appendix  C.  It  is  also  necessary  to  be  ab'e  to  numerically 
convolve  probability  distributions  in  this  analysis.  The  Fast 
Fourier  Transform  is  used  to  do  this.  This  procedure  is 
described  in  Appendix  D.  As  part  of  the  probability  of 

horizontal  overlap  estimation  one  needs  to  estimate  the 

probability  of  alongtrack  separation.  This  estimate  is 

developed  in  Appendix  E. 

In  all  of  the  analyses  mentioned  above  one  needs  the  probability 
that  an  aircraft  pair  has  a  particular  crosstrack  separation  and 
crosstrack  closing  speed.  Since  the  data  is  taken  on  single 

aircraft,  a  procedure  is  required  to  convert  the  single  aircraft 
data  into  aircraft  pair  data.  This  procedure  is  described  in 
Appendix  F. 

The  controller  intervention  rate  due  to  conflict  alerts  was 
estimated  by  simulation.  The  simulation  is  described  in 

Appendix  G.  The  NAS  Conflict  Alert  function  which  was  emulated 
in  the  simulation  is  described  in  Appendix  H. 

A  list  of  all  the  variables  and  symbols  used  in  both  this  volume 
and  Volume  I  can  be  found  in  the  glossary  in  Appendix  I.  The 
references  for  this  volume  are  found  in  Appendix  J. 
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APPENDIX  A 

THE  CONFLICT  REGION  BOUNDARY 


This  appendix  will  develop  the  concept  o£  the  conflict  region 
boundary.  As  indicated  in  Section  3  of  Volume  I,  the  conflict 
region  boundary  in  the  analysis  is  a  straight  line  on  the 
crosstrack  separation,  crosstrack  closing  speed  plane.  The 
equation  for  this  line  and  the  assumptions  that  were  made  in 
constructing  the  line  will  be  developed  in  this  appendix. 

The  situation  is  the  following:  A  pair  of  aircraft  are  flying 
the  same  direction  on  two  parallel  routes.  They  are  assumed  to 
have  the  same  forward  speed  but  each  aircraft  can  have  different 
alongtrack  and  crosstrack  component  speeds.  The  aircraft  are 
positioned  on  their  respective  routes  as  shown  in  Figure  A-l. 
The  aircraft  are  projected  ahead  along  straight  line  paths  for  a 
time  TL.  If  within  that  time  they  are  separated  by  less  than  a 
distance  D,  then  the  aircraft  pair  is  said  to  be  in  potential 
conflict. 

The  sign  conventions  for  the  deviations  and  velocities  are  as 
follows:  Crosstrack  deviations  and  crosstrack  velocities  for 

the  individual  aircraft  are  measured  with  respect  to  the  route 
centerlines  of  the  respective  routes  and  are  positive  in  the 
direction  toward  the  other  route.  In  Figure  A-l,  and  Vy£ 

are  the  crosstrack  deviation  and  velocity,  respectively,  for  an 
aircraft  on  route  L.  The  crosstrack  separation  between  a  pair 
of  aircraft  is  denoted  as  y  and  is  positive  when  the  sum  of  the 
crosstrack  deviations  is  less  than  the  spacing  between  the 
routes.  Alongtrack  displacement  between  the  aircraft,  x,  is 
measured  with  respect  to  the  aircraft  on  route  1.  The 
alongtrack  displacement  x,  is  positive  when  the  aircraft  on 
route  2  is  ahead  of  the  aircraft  on  route  l.  With  these 
conventions,  the  distance  between  the  aircraft  (d)  as  a  function 
of  time  (t)  (assuming  rectilinear  motion)  is 

“  -  \f  +  (v*2  -  v) '} 2  +  {>  -  (v +  v) c}2  (A-i) 

If  we  let 


x  *  VX2  *  vxl  (alongtrack  closing  speed) 

• 

y  y2  +  V  (crosstrack  closing  speed) 


A-l 


and  solve  (A-l)  for  y  then 


fl 


y  -  yt  ±\[d2  -  x2  -  t 2  x2  -  2  txx 


(A-2) 


The  objective  in  this  analysis  is  to  arrive  at  a  relationship 
between  the  crosstrack  separation  of  an  aircraft  pair  and  the 
conditions  under  which  a  potential  conflict  would  exist.  Thus, 
one  would  want  to  know  when  the  aircraft  pair  first  enters  into 
potential  conflict.  This  can  be  determined  by  letting  d  equal 
the  threshold  distance  for  a  potential  conflict,  D,  and  t  equal 
the  look-ahead  time,  TL.  Hence,  one  finds  the  crosstrack 
separation,  y,  at  which  the  threshold  separation,  D,  will  be 
achieved  at  time  TL  in  the  future. 

It  is  assumed  at  this  point  that  prior  to  entering  into 
potential  conflict  the  aircraft  are  oriented  with  respect  to 
each  other  in  the  same  way  that  their  respective  routes  are 
oriented  and  they  are  closing  in  the  crosstrack  direction.  In 
other  words,  with  reference  to  the  sign  convention  used  in 
Figure  A-l,  the  crosstrack  separation  and  the  crosstrack  closing 
speed  are  both  positive  prior  to  entering  into  potential 
conflict.  With  this  assumption,  the  plus  sign  in  front  of  the 
radical  in  equation  A-2  is  chosen. 

Making  the  substitutions  into  equation  A-2  we  can  define  the 
conflict  region  boundary  as: 


y 


2 

x 


.2 

x 


-  2TL  x 


x 


(A-3) 


Given  TL,  D,  x,  and  x  one  could  plot  the  conflict  region 
boundary  on  the  y  plane.  However,  x  in  general  can  take  on 
many  values  for  a  given  value  of  y.  In  fact,  if  the  maximum 
crosstrack  speed  of  an  aircraft  is  then  the  range  of 

alongtrack  closing  speeds  (&)  for  a  given  crosstrack  closing 
speed  ($)  is 


ymax 


(w  -V^ 


2 

ymax 


T 


where  V  is  the  forward  speed  of  the  two  aircraft.  It  should  be 
remembered  that  the  x  dimension  is  measured  with  respect  to  the 
position  of  the  aircraft  on  route  l. 

If  it  is  assumed  that  the  distributions  of  crosstrack  speed  for 
single  aircraft  on  the  two  routes  are  the  same  and  the  aircraft 
all  have  equal  speed,  then  the  expected  alongtrack  closing  speed 
(R(fc))  over  all  values  of  y  will  be  zero.  To  get  a  first  order 
look  at  the  conflict  boundary  we  substitute  x=0  into  (A-3)  to 
arrive  at 


y 


y  TL  +\£2" -  x2 


(A-5) 


The  conflict  region  based  on  (A-5)  is  shown  in  Figure  A-2. 
Recall  that  the  line  in  Figure  A-2  defines  the  point  at  which 
the  threshold  separation,  D,  is  projected  to  be  violated  at  a 
time  TL  in  the  future.  Since  the  aircraft  pair  is  closer 
together  for  smaller  values  of  y,  the  region  to  the  left  of  the 
line  in  Figure  A-2  represents  combinations  of  y  and  y  where  the 
threshold  separation  is  violated  prior  to  time  TL.  Note  also 
that  the  y-intercept  of  this  curve  is  ^D2  I  x2^  which 
depends  on  the  initial  alongtrack  separation  x.  At  the  point 
y=0  the  horizontal  separation  will  be  ^/D2  -  x2.  In  other 

words,  if  the  aircraft  were  not  closing  on  each  other  (v=0.x=0) 
but  their  horizontal  separation  were  less  then^/D2  -  x2  then 
the  aircraft  pair  would  be  in  potential  conflict.  The  maximum 
value  for  y  is  2  Vy^a„  since  each  aircraft  could  be  flying 
away  from  its  route  centerline  with  no  more  than  the  maximum 
crosstrack  speed  Vy^*. 

When  proximity  is  discussed  in  the  Conflict  Monitoring  Analysis, 
its  meaning  is  based  on  the  allowable  range  of  x  values.  One 
thing  which  should  be  noted  about  equation  (A-5)  is  that  there 
is  a  limited  range  for  x.  If  |x|>D,  then  the  square  root  in 
equation  (A-5)  gives  an  imaginary  number.  This  means  that  if 
x=0  and  the  aircraft  are  separated  alongtrack  by  more  than  D 
then  there  can  be  no  potential  conflict.  In  other  words,  if  the 
aircraft  are  not  closing  in  the  alongtrack  direction  (x**0),  then 
no  matter  how  close  they  are  in  the  crosstrack  direction  they 
will  always  be  separated  by  more  than  the  distance  D.  Aircraft 
which  are  spaced  alongtrack  greater  than  a  distrance  D  are  not 
proximate  because  they  cannot  be  in  potential  conflict  and  hence 
the  analysis  indicates  that  they  cannot  collide. 
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THE  CONFLICT  REGION 


The  expression  for  the  conflict  region  boundary  given  in 
equation  (A-S)  and  illustrated  in  Figure  A-2  was  used  in  the 
Conflict  Monitoring  Analysis.  Recall  that  it  was  derived  from 
equation  (A-3)  by  replacing  x  by  its  expected  value.  Equation 
(A-S)  was  used  because  it  is  simpler  than  (A-3)  (being  a 
straight  line)  and  it  is  a  reasonable  approximation  to  the 
conflict  region  boundaries  with  x+0 . 

To  see  how  good  an  approximation  equation  (A-5)  is  we  shall  go 
back  and  investigate  equation  (A-3).  From  equation  (A-3)  it  is 
apparent  that  for  a  given  value  of  y,  the  value  of  y  will  be  a 
minimum  when 


z 


(X  +  Tl  i)2 


(A-6) 


is  a  minimum  greater  than  or  equal  to  zero.  The  value  of  y  will 
be  a  maximum  when  z  is  a  maximum.  However,  it  is  necessary  from 
(A-4)  that 


<  x<'x2 


(A-7) 


where 


and 


(y 


*> 

V  )‘ 

ymax 


If  we  let  M  ■  x  ♦  TL  x,  then  equation  (A-6)  is  a  parabola  as 
shown  in  Figure  A-3.  The  value  of  M  must  be  between  the  values 
of  and  M2  to  satisfy  equation  (A-2).  If  either  M}  or 

K»  is  outside  the  range  -D£M£D  then  the  minimum  s  is  zero. 
If  both  Mj  and  M2  are  iqu.de  the  range  -D$  H£  D  then  z  is 
the  minimta  of  D2  -  and  D2  -  M?.  T3ie  maximum  z 

is  the  maximum  of  D2  -  Mj  and  D2  “-  M(  if  Mj  and 

M2  are  both  on  the  same  side  of  the  origin,  "if  M}  and  M2 
are  on  opposites  sides  of  the  origin  the  maximia  z  is  D2, 
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MAX/MIN  REGION  FOR 
DISCRIMINANT  IN  EQUATION  A-3 
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From  the  above  considerations  the  conflict  region  boundaries  for 
the  range  of  possible  x  values  given  in  equation  (A-7)  can  be 
drawn  for  various  values  of  alongtrack  separation  (x).  Figures 
A-4  through  A-6  show  a  set  of  conflict  region  boundary  envelopes 
for  x  *  0.0  nmi,  x  “  +3.0  nmi,  and  x  =  +5.0  nmi.  If  a  truly 
conservative  risk  estimate  were  to  be  made,  then  the  left-most 
conflict  region  boundary  should  be  used.  This  boundary  would 
allow  the  aircraft  pair  to  be  closer  together  before  the 
resolution  maneuver  is  executed  and  thus  increase  the  risk  of 
collision. 

To  get  an  appreciation  of  the  quality  of  the  use  of  x  =  0  in 
equation  (A-3)  it  is  instructive  to  look  at  the  distributions  of 
early  and  late  detections  which  would  result  from  the  use  of  the 
x  *  0  conflict  region  boundary.  Such  distributions  were 
constructed  for  a  number  of  values  of  the  crosstrack  closing 
speed  (y)  for  several  values  of  x.  A  set  of  these  distributions 

are  shown  in  Figures  A-7,  A-8,  and  A-9.  Each  distribution  is 

scaled  relative  to  the  frequency  with  which  the  particular  value 
of  y  is  expected  to  be  observed. 

Figure  A-7  shows  that  when  x  *  0,  the  aircraft  are  nearly  always 
detected  early  using  the  x  =  0  conflict  region  boundary. 

However,  the  maximum  time  difference  is  only  a  few  seconds.  In 
Figure  A-8,  when  x  =*  +3  nmi,  the  aircraft  pairs  could  be 

detected  early  or  late.  At  the  lesser  values  of  y  the 

distribution  is  more  or  less  symmetrical,  while  for  greater 
values  of  y  the  distribution  is  skewed  toward  a  later  (more 
conservative)  detection.  When  x  =  +5 nmi  (see  Figure  A-9),  one 
always  detects  late  by  using  the  x  =  0  conflict  region 

boundary.  The  spike  at  y  5  10  kts  indicates  that  20X  of  the 

aircraft  pairs  will  be  detected  more  than  200  seconds  late. 
Thus,  from  the  consideration  of  the  early/late  detection 

distributions  it  appears  that  the  use  of  the  x  ■  0  conflict 
region  boundary  is  prudent  if  not  some  what  conservative. 

Another  consequence  of  not  letting  x  *  0  is  that  x  would  no 
longer  be  required  to  be  within  the  range  -Dsx  SD  to  define  a 
conflict  region  boundary.  If  x  i  0  and  |x|  is  chosen  greater 
than  D,  there  will  be  a  range  of  conflict  region  boundaries  as 
dep  ted  in  Figure  A- 10.  This  figure  represents  the  envelope  of 
conflict  region  boundaries  for  |x|  =  6.5  nmi.  For  |x|  greater 

than  this  the  aircraft  would  be  too  far  apart  to  be  in  conflict 
even  when  x  ^  0. 
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ENVELOPE  OF  CONFLICT  REGION  BOUNDARIES 
FOR  x  =  ±5.0  NMI 
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DISTRIBUTION  OF  EARLY  (-)  OR  LATE  (  +  )  DETECTIONS 
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FOR  x  =  ±6.5  NMI 


One  assumption  made  in  this  analysis  should  be  highlighted.  For 
this  analysis  the  forward  speeds  of  both  aircraft  were  assumed 
to  be  the  same.  For  the  high  altitude  CONUS  airspace  with  many 
aircraft  with  the  same  performance  capabilities  flying  at  the 
same  altitude  this  assumption  should  be  a  reasonable.  However, 
there  are  some  slower  aircraft  in  the  high  altitude  airspace. 
If  one  wanted  to  account  for  these  speed  differences,  the  entire 
analysis  could  conceivably  be  done  on  a  pair-by-pair  basis  with 
the  results  properly  weighted.  For  expediency,  the  different 
speed  case  was  not  considered  in  this  analysis  but  will  be 
addressed  in  the  future. 
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APPENDIX  _B 

COMPUTATION  OF  THE  PROBABILITY  OF  HORIZONTAL  .OVERLAP 
Bl.  INTRODUCTION 


In  order  for  a  pair  of  aircraft  to  collide,  the  pair  must  first 
penetrate  the  conflict  region  boundary  described  in  Appendix  A. 
Therefore,  the  probability  of  horizontal  overlap  can  be 
calculated  as  the  product  of  the  probability  that  the  aircraft 
pair  is  on  the  conflict  region  boundary  and  the  conditional 
probability  that  the  pair  has  a  delay  time  and  an  avoidance 
maneuver  turn  rate  which  would  result  in  horizontal  overlap 
given  that  the  pair  is  on  the  conflict  region  boundary. 

First,  this  appendix  will  discuss  the  distributions  of  the  time 
delay  and  the  turn  rate  which  were  used  lor  this  analysis. 
Second,  the  appendix  will  describe  using  these  distributions  to 
calculate  the  conditional  probability  that  an  aircraft  pair  has 
a  delay  time  and  a  turn  rate  which  would  result  in  horizontal 
overlap  given  that  the  pair  is  on  the  conflict  >egion  boundary. 
Finally  this  appendix  will  illustrate  the  calculation  of  the 
probability  that  an  aircratt  pair  is  on  "he  conflict  region 
boundary.  This  then  provides  all  the  informal'  in  necessary  to 
estimate  the  probability  of  horizontal  overlap. 

B2.  DELAY  TIME 


The  delay  time  in  the  model  is  defined  to  he  that  time  when  both 
aircraft  are  in  the  conflict  region  and  ate  flying  straight. 
This  includes  the  time  between  the  actua'  penetration  of  the 

conflict  region  boundary  and  the  detection  of  the  potential 
conflict.  The  delay  time  will  also  include  the  time  taken  for 
the  controller  to  recognize  the  conflict,  to  decide,  and  to 

communicate  instructions  to  the  pilot,  and  for  the  pilot  to 
start  to  take  action  on  the  instructions.  The  total  delay  time 
will  be  the  sum  of  the  above  listed  delays.  We  will  first 
examine  the  delay  due  to  the  detection  process  and  then  address 
the  remainder  of  the  delay. 

B2.1  Delay  Due  to  Detection 

Consider  a  pair  o  aircraft  transgressing  the  conflict  region 
boundary  as  shown  in  Figure  B-l.  At  radar  update  1  the  pair  is 
just  outside  the  conflict  region  boundary.  At  update  2  the  pair 
is  just  inside  the  boundary.  At  subsequent  updates  the  pair  is 

further  inside  the  conflict  region.  As  shown  in  the  figure,  it 
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FIGURE  B-1 

PROGRESSION  OF  A  PAIR  OF  AIRCRAFT 
ACROSS  THE  CONFLICT  REGION  BOUNDARY 


is  assumed  Chat  Che  crosstrack  closing  speed  is  constant.  This 
assumption  is  based  on  the  prior  assumption  that  the  velocities 
of  the  individual  aircraft  remain  constant  after  the  pair 
crosses  the  conflict  region  boundary.  As  discussed  in  Appendix 
A  the  conflict  region  boundary  is  defined  under  the  assumption 
that  x»0.  In  terms  of  Figure  B-l  this  assumption  means  that  on 
the  average  the  conflict  region  boundary  is  the  same  for  each 
radar  update.  Even  though  the  expected  value  of  x  is  zero,  it 
can  have  values  in  the  range  given  in  equation  (A-4).  For  most 
of  the  values  of  4  in  this  range  the  change  in  the  value  of  x 
over  a  update  interval  of  time  will  be  insignificant.  The  cases 
of  larger  changes  of  x  over  an  update  interval  will  be  very 
infrequent.  For  a  given  value  of  x  (alongtrack  separation)  and 
y  (crosstrack  closing  speed)  we  want  to  construct  a  probabi  lity 
density  function  of  the  time  interval  between  the  time  the 
aircraft  pair  enters  into  the  conflict  region  and  the  time  it  is 
detected  as  being  in  potential  conflict. 

The  detection  process  is  based  on  a  set  of  surveillance  returns 
which  have  been  processed  through  an  alpha-beta  tracker.  The 
position  and  velocity  estimates  from  the  tracker  are  assumed  to 
have  normally  distributed  correlated  errors  with  zero  bias. 
Referring  to  Figure  B-2,  one  can  visualize  the  detection 
process.  The  aircraft  pair  crosses  into  the  conflict  region 
with  a  given  crosstrack  closing  speed  (y)  and  is  a  distance  F 
inside  the  conflict  region  boundary  when  the  initial  radar 
observation  is  made.  After  the  next  revolution  of  the  radar 
antenna  (after  dt  hoursl  ,  the  aircraft  pair  has  a  crosstrack 
separation  of  Y2  F*  The  crosstrack  separation  lost  between 
Y^p,  and  Y2p  is  G-ydt. 

At  each  crosstrack  separation  Yj  p,  there  is  a  probability 
PNDj  p  that  the  estimates  of  the  crosstrack  separation  and 
closing  speed  will  indicate  that  the  aircraft  pair  is  still 
outside  the  conflict  region.  The  computation  of  this 
probability  is  discussed  in  Appendix  C. 

Since  each  radar  update  is  assumed  to  give  an  independent 
measurement  of  the  aircraft  positions,  the  probability  of  not 
being  detected  in  the  conflict  region  during  the  first  i 
updates  is  given  as 
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FIGURE  B-2 

CROSSTRACK  SEPARATIONS  DURING  THE  DETECTION  PROCESS 
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The  probability  of  first  detecting  an  aircraft  pair  in  the 
conflict  region  in  the  class  interval  i  (see  Figure  B-2 1  is 
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PFDi,F  *  PND  i_l,F  (l-PNDifF)  (B-?) 

It  should  be  remembered  that  the  initial  observation  of  aircraft 
pair  was  made  when  the  aircraft  pair  was  inside  the  conflict 
region  by  a  crosstrack  separation  distance  F.  The  distance  F  is 
a  random  variable  with  a  probability  density  function  q(F). 
Since  the  time  of  entering  the  conflict  region  is  in  no  way 
related  to  the  timing  of  the  radar  scans,  it  is  reasonable  to 
assume  the  F  is  uniformly  distributed  between  0  and  G.  Thus 

q(F)=l/G  0  <  F  <  G  (B-3) 


Therefore,  the  probabilty  of  first  detection  in  class  interval  i 
is 
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In  the  actual  computation  a  set  of  N  equally  spaced  initial 
positions  were  selected.  Thus 

n 

PFDi  3  NG~  ^  PFDi,(j-l)G/N  (B-3) 

The  quantity  computed  in  (B-5)  is  the  probability  that  an 
aircraft  pair  with  crosstrack  closing  speed  y  will  be  first 
detected  in  class  interval  i.  The  class  interval  width  is  G 
(see  Figure  B-2)  is  the  crosstrack  separation  lost  during  one 
radar  update. 

2  Delay  Due  to  the  Controller,  the  Communication  Link,  and  the 
Pilot 


After  the  aircraft  pair  is  detected  as  being  within  the  conflict 
region  there  will  be  an  additional  delay  while  the  controller 
decides  what  to  do  and  the  resolution  commands  are  transmitted 
to  the  pilot.  The  pilot  then  has  to  decide  what  to  do  and  start 
to  turn  his  aircraft. 
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The  various  components  of  this  delay  are  very  difficult  to  model 
because  of  the  human  element.  However,  there  is  data  available 
from  a  simulation  done  in  Great  Britain  on  controller  response 
to  threshold  transgression.1^  The  histogram  of  delay  times 

with  a  computer-assisted  system  is  shown  in  Figure  B-3.  These 
delay  times  were  measured  from  the  time  of  presentation  to  the 
controller  to  the  time  the  controller  communicates  with  the 
pilot. 

Even  though  the  controllers  in  the  simulation  had  a  workload  in 
addition  to  their  monitoring  role,  one  must  still  remember  that 
the  controllers  knew  that  they  were  in  a  simulated  environment. 
For  this  reason  the  histogram  data  in  Figure  B-3  were  fit  to  a 
Ganna  function  to  give  the  delay  distribution  a  long  tail.  The 
Gamma  function  fit  is  also  shown  in  Figure  B-3. 

The  collision  risk  will  be  sensitive  to  the  delay  time.  This 

means  that  the  conservativeness  of  the  analysis  could  be 
dictated  by  the  chosen  delay  function.  The  particular  delay 

function  which  was  fit  to  the  histogram  data  has  a  long  tail 
which  was  truncated  at  600  seconds  for  computational  reasons. 
This  length  of  time  is  approximately  one-half  the  flying  time 
through  a  sector.  With  this  length  of  delay  one  could  argue 

that  possible  failure  of  the  controller  to  notice  the  potential 
conflict  or  short  term  communication  breakdowns  or  other  outages 
are  essentially  accounted  for. 

B2.3  Total  Delay 

The  delay  due  to  the  detection  process  is  independent  of  the 
delay  due  to  the  controller/communication/pilot  reactions. 
Therefore,  the  total  delay  is  the  sum  of  two  independent  random 
variables.  The  probability  density  function  of  the  sun  of  two 
independent  random  variables  is  the  convolution  of  their 
probability  densities.  The  technique  used  to  convolve  the 
probability  density  functions  is  the  Fast  Finite  Fourier 
Transform.  The  Fast  Finite  Fourier  Transform  technique  is 
discussed  in  detail  in  Appendix  D. 

B.3-  .TffiB-RAIE 

In  the  Conflict  Monitoring  Analysis  the  turn  rate  is  specified 
as  a  range  of  bank  angles  since  it  is  the  bank  angle  that  the 
pilot  controls  when  he  makes  a  turn. 

Wte  will  assune  that  the  bank  angle  chosen  by  the  pilot  to  make 
his  turn  back  to  his  assigned  route  centerline  will  be  in  the 
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FIGURE  B-3 
DELAY  DISTRIBUTION 
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range  of  10  degrees  to  10  degrees.  The  minimum  Hank,  angle  of  10 
degrees  means  that  the  pilot  will  make  a  definite  resolution 
maneuver.  The  maximum  bank  angle  of  10  degrees  was  chosen 
because  it  is  the  angle  which  usually  demonstrates  the  threshold 
of  passenger  discomfort  due  to  g  forces.  It  is  also  assumed 
that  the  choice  of  each  bank  angle  within  this  range  is  equally 
likely.  The  probability  density  function  of  the  bank  angle  is 
shown  in  Figure  B-4. 

The  distribution  of  bsnk  angles  must  be  related  to  a 
distribution  of  turn  rates  in  order  to  find  the  relationship 
between  the  turn  and  the  time  and  distance  required  for  the 
maneuver.  If  the  pilot  makes  a  coordinated  turn  then 
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where  uj  is  the  turn  rate,  g  is  the  gravitational  constant,  V  is 
the  forward  velocity  of  the  aircraft,  and  *  is  the  bank  angle. 


If  the  pdf  of  k  is  uniform  as  shown  in  Figure  B-4,  then  the  pdf 
of  can  be  shown  to  be 
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The  pdf  j(w)  that  corresponds  to  the  bank  angle  pdf  in  Figure 
B-4  is  shown  in  Figure  B-5. 


B4.  HORIZONTAL  OVERLAP  REGION 


Up  to  this  point  we  have  developed  the  probability  density 
function  of  the  total  delay  after  an  aircraft  pair  penetrates 
the  conflict  region  and  the  probability  density  function  of  the 
turn  rate  used  by  an  aircraft  to  return  toward  its  route  center 
line.  The  problem  now  reduces  to  finding  the  probability  that 
the  aircraft  pair  will  come  into  horizontal  overlap  given  that 
1)  the  pair  passes  into  the  conflict  region,  2)  there  is  a 
delay,  and  3)  there  is  a  turn  back  toward  the  route  centerline 
by  one  of  the  aircraft. 

After  the  aircraft  pair  passes  into  the  conflict  region,  the 
aircraft  are  assumed  to  continue  flying  in  a  straight  line  for  a 
delay  time  t^.  Then  one  aircraft  will  make  a  horizontal  turn 


with  a  turn  rate  w  •  At  some  point  after  the  pair  enters  the 
conflict  region,  either  a  minimum  separation  between  the 
aircraft  will  he  achieved  after  which  the  aircraft  will 

separate,  or  the  aircraft  will  collide.  With  a  delay  and  then  a 
turn  there  are  four  ways  in  which  minimum  separation  can  occur. 
First,  the  minimum  separation  can  occur  at  that  time  when  the 

velocities  of  the  two  aircraft  become  parallel.  A  second  way  is 
for  the  aircraft  to  reach  their  minimum  separation  before  one  of 
the  aircraft  starts  to  turn.  A  third  way  would  be  for  the 

minimum  separation  to  occur  during  the  turn  and  before  the 
velocities  become  parallel.  The  fourth  case  is  where  the 

absolute  minimum  separation  occurs  after  the  velocities  become 
parallel.  Only  the  first  three  cases  are  of  interest  because 
the  fourth  case  implies  that  the  paths  of  the  two  aircraft  cross 
without  a  collision  and  then  one  of  the  aircraft  turns  back  into 
the  path  of  the  other  aircraft. 

To  analyze  the  first  three  cases  we  will  start  the  aircraft  pair 
in  the  condition  of  horizontal  overlap.  The  aircraft  will  then 
be  "flown  backwards  in  time".  One  aircraft  will  execute  a  turn 
at  a  turn  rate  u  while  the  other  flies  straight.  The  turning 
aircraft  will  then  come  out  of  its  turn  and  fly  a  straight 
path.  Both  aircraft  fly  (backwards)  along  their  respective 
straight  paths  until  they  reach  positions  which  will  place  them 
on  the  boundary  of  the  conflict  region.  By  doing  the  analysis 
this  way  one  can  identify  those  combinations  of  turn  rate  and 
delay  time  which  would  place  an  aircraft  pair  A;n  horizontal 
overlap  given  that  the  pair  started  from  a  position  on  the 
conflict  region  boundary. 

To  find  these  combinations  of  delay  times  and  turn  rates  we 
consider  the  scenario  shewn  in  Figure  B-b.  Aircraft  1  is  on  the 
conflict  region  boundary  at  position  Xjq  Yjq  with  velocity 
vxl,  vyl  while  aircraft  2  is  at  position  *20,  Y20  with 
velocity  Vx2,  Vy2- 

The  velocities  are  such  that 
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The  trajectory  of  aircraft  2  which  flies  the  straight  course  is: 
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The  trajectory  of  aircraft  1  which  flies  the  straight  course  for 
a  time  tj  then  a  curved  course  for  a  time  tt  is: 
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The  time  that  aircraft  1  is  in  the  turn,  tc,  ir  given  by 

*  . 


t  =■  p  -  -  ■  -  ;  o<p<i 

C  C J 


(B-Jl) 


where  P  is  a  factor  which  delineates  the  three  cases  discussed 
above.  If  P“0  then  there  is  no  turn  and  the  minimum  separation 
will  be  when  both  aircraft  are  flying  straight.  If  P=1  then 
aircraft  1  will  have  made  its  turn  to  come  parallel  to  aircraft 


2  at  the  point  of  minimum  separation.  If  0  i  P  •  1  then  the 
separation  at  time  t  will  be  during  the  turn  of  aircraft  1 
before  the  velocities  of  the  two  aircraft  become  parallel. 

There  are  two  conditions  that  are  imposed  on  this  scenario: 

1.  The  aircraft  pair  starts  on  the  conflict  region 
boundary 

2.  The  aircraft  pair  ends  up  in  a  horizontal  overlap 
situation. 

The  first  condition  is  specified  by  the  expression 


Y20'Y10  "  <Vyl-Vy2)TL  ♦  VPhX20  - 


where  TL  is  the  look-ahead  time  and  D  is  the  threshold 
separation  which  is  tested  for  in  the  conflict  alert  function. 
The  second  condition  is  that  the  aircraft  pair  end  up  in 
horizontal  overlap.  If  we  assume  that  the  aircraft  are 
represented  as  right  cylinders  with  radius  R,  this  means  that 
the  aircraft  will  be  in  overlap  if  their  centers  are  separated 
by  less  than  a  distance  2R.  Therefore,  if 


jxjCtl-Xj/tlj.  +  -|y 2  ( 1 1  -  Y T ( t  )J-  e  (2R)2 


(B-l.l) 


then  the  aircraft  will  be  in  a  conditon  of  horizontal  overlap  at 
time  t. 

Substituting  (B-9)  and  (B-10)  into  (B-13)  and  using  (B-8), 
(B-ll),  and  ( B— 12)  we  can  rearrange  terms  to  arrive  at 


td2(I2  +  r2)  +  td  ( 2Ix  +  2IJH  +  2KLH  +  2  KM)  +  (x2  + 

J2  H2  +  2 JHx  +  L2  H2  +  M2  +  2LMH  -  4R2)-£0  (B-14) 
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vherp  H3!/-, 

I  ™  Vf  a i n  v  -  sinf) 

J  *  V  [P(v-£)sin  i’  -  cns  (.  *■  cost  4  ( 1  ~ P )  +  Pr )] 

K  3  V  ( cos  v  -  cob  £,  ) 

L  3  v[P(  v-£)  cos  v  sin(_4(l-P)  +  Pv)  +  sin4] 

M  -  fVyl  -  Vy2)  TL  +  v/D2  -  x7  =  y 

One  can  find  the  range  of  values  of  tj  for  which  equation 
(B-14)  ia  true  by  using  the  quadratic  formula.  Equation  (B-14) 
will  be  true  for  tj  between  the  roots  of  the  expresion 
obtained  by  setting  the  lefthand  side  of  equation  (B-14)  equal 
to  zero.  However,  in  order  that  this  expression  have  real  roots 
the  following  expression  must  be  true: 


-H2(il  -  JK)2  +  2H(IL  -  KJ)  (Kx  -  IM)  +  4R2  (I2  +  K2) 
+  2 IKMx  -  I2M2  -  K2  x2  >0 


Therefore,  restricting  our  attention  to  the  range  of  values  of  H 
for  which  (B-15)  is  true  wiLl  insure  that  there  will  be  real 
values  of  t<j  which  satisfy  (B-14).  However,  in  order  that 
(B-15)  be  true  for  real  values  of  H,  H  must  lie  between  the 
roots  of  the  expression  obtained  by  setting  the  lefthand  side  of 
(B-15)  equal  to  zero.  This  expression  will  have  real  roots  if 
the  following  inequality  holds: 


(IL-KJ) 2 (Kx-IM) 2  +  (IL-JK)2(4R2(I2+k2)  +  2 IKMx 


9  9  9  9 

I'-'T-irx-) 


>  0 


But  (B-16)  reduces  to 


(B-16) 


R2  (IL  -  KJ)2  (I2  +  K2)>0  (B-17) 

which  is  always  true.  Therefore,  the  solution  to  (B-15)  is  for 
H  between 


(Kx-IM)  -  2R 
IL  -  JK 


B- 1  r> 


and 


(Kx  -  IM)  +  2R 
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+  K 


IL  -  JK 


The  additional  stipulation  on  H  is  that  it  be  positive.  This 
means  that  the  turning  aircraft  is  turning  back  towards  its 
assigned  route  rather  than  away  from  it.  This  is  the  range  of  H 
values  for  which  (B-14)  is  satisfied  by  real  values  of  t<j.  As 
an  example  of  three  different  conditions  consider  Figure  B-7 . 
Given  the  radius  of  the  aircraft  (R),  the  initial  alongtrack 
spacing  (x),  and  the  initial  crosstrack  closing  speed  (y),  a 
"crescent-shaped"  region  in  t^- w  space  can  be  constructed. 
The  region  inside  the  crescent  are  those  t,j,u>  combinations 
(given  R,x,  and  y)  for  which  the  aircraft  pair  start  on  the 
conflict  region  boundary  and  end  in  a  horizontal  overlap 
c  ond  i  t  i  on . 

To  compute  the  probability  of  being  within  a  crescent  shaped 
region  in  Figure  B-7  one  would  really  have  to  ask  for  the 
probability  of  also  being  within  Ax  of  x  and  Ay  of  y.  To 
approach  this  problem  numerically  we  will  divide  up  the  range  of 
y  and  the  range  of  x  into  small  cells.  For  instance,  consider 
Figure  B-8.  Here  we  have  a  cell  which  is  1  nmi  in  x  and  10 
knots  in  y.  If  we  were  to  select  those  x,  y  combinations  shown 
by  the  x's  in  Figure  B-8,  we  could  construct  a  t<j-w  crescent 
region  for  each.  The  result  of  doing  this  for  P=1  and 

Vy  1-200  kts  is  shown  in  Figure  B-9.  On  the  scale  shown  in 
Figure  B-9,  the  crescent  regions  are  lines. 

If  other  values  of  P  between  0  and  1  and  other  values  of 
are  used,  Figure  B-9  is  expanded  to  get  a  region  in  t^-w  space 
such  as  the  one  shown  in  Figure  B-10.  This  figure  was 
calculated  for  P*0,l/2,  and  1,  and  for  five  different  values  of 
Vyi  and  Vy2  which  could  result  in  y.  The  dots  represent  the 

length  and  width  of  each  crescent  shaped  region. 

The  dots  in  the  upper  left  hand  corner  of  Figure  B-10  are  those 
which  correspond  to  the  minimum  time  delay  necessary  for  the 
centers  of  the  aircraft  to  overlap  when  no  turn  is  executed. 

One  way  in  which  these  points  can  be  identified  is  by  performing 
a  grid  search  over  various  values  of  crosstrack  closing  speed 
(y),  alongtrack  separation  (x),  and  crosstrack  speed  partition 

between  the  two  aircraft  (F,  to  be  defined  below). 
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It  could  be  that  a  particular  x,  y  cell  does  not  have  any 
"atraight  line"  overlaps.  The  problem  is  to  determine  the 
existence  of  such  straight  line  overlaps  within  an  x,  y  cell  and 
if  there  are  such  points  then  find  the  one  with  the  minimum  time 
delay  from  the  trangression  of  the  conflict  region  boundary  to 
overlap. 

This  problem  is  a  nonlinear  programming  problem  with  equality 
and  inequality  constraints.  The  problem  can  be  structured  as: 


I 


minimize  tj  ■  TL  +  v/d2  -  X2 


subject  to: 


F-l  <0 


-  1-F  £  0 


y  -  yi<  o 
Y2  -  y so 
x  -  xx  <o 

X2  -  X  <  o 

( l+B2)x2  -  2TLyBx  +  TL2y2B2  -  D2B2  =  0 


(B-18) 


(B-I9) 


where 


/V2  -  (FVm-y)2 


TL  *  look  ahead  time 
D  =  Minimum  radar  separation 

Vm  =»  Maximum  observed  crosstrack  velocity  (single 
aircraft) 

xl»x2  *  limits  on  x 
^1>Y2  *  limits  on  $ 


The  object  is  to  find  x,  y,  and  F  such  that  the  above  conditions 
are  satisfied.  Since  the  general  solution  is  difficult  we  will 
investigate  the  existence  of  straight  line  overlap*  by 
considering  equation  (8-19).  This  equation  sets  the  condition 
that  the  aircraft  overlap.  Solving  this  equation  for  x  gives: 
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/  2 
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x  -  Ti.  y  B  ♦  1 

V  D  (1 

♦  B  )  -  B 
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1 

+  B 
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If  x  is  substituted  back,  into  the  expression  from  which  equation 
(B-19)  was  derived,  it  i9  found  that  the  only  applicable  root 
for  x  is  the  one  with  the  sign  in  equation  (B-10).  The 

values  of  x  for  the  spectrum  of  F  and  v  values  are  shown  in 
Figure  B-ll.  From  this  figure  one  can  see  that  for  cells  with 
I X |  >3.5  nrai  there  are  no  straight  line  overlaps  --  aircraft  are 
too  far  apart  from  each  other  in  the  alongtrack  direction  for 
there  to  be  an  overlap.  Once  the  existence  of  a  straight  line 
overlap  is  determined  for  a  particular  x,  v  cell,  the  problem  is 
to  find  the  combination  of  values  of  x,  y,  and  F  which  minimize 
fcd  given  in  equation  (B-19).  A  close  approximat ion  to  the 
minimum  value  of  t<j  is  to  choose  the  maximum  value  of  y  in  the 
cell  and  then  choose  the  maximum  feasible  value  of  x  given  that 
y.  This  then  is  how  the  point  labeled  "A"  in  Figure  B-10  is 
determined. 

For  any  time  delay  greater  than  A  for  this  x,  y  cell,  the 
assumption  is  made  that  the  aircraft  pair  will  overlap  before 
starting  an  avoidance  maneuver.  For  time  delays  from  A  down  to 
where  the  pattern  of  dots  ends,  one  would  say  that  there  was 
horizontal  overlap  for  those  time  delays  at  the  particular  turn 
rate.  To  describe  this  region  an  envelope  of  the  lower  boundary 
of  the  pattern  of  dots  was  approximated  by  a  function  of  the 
form: 


t j*A  -  C/w. 


(B-21) 


The  curved  dashed  lines  in  Figure  B-10  show  the  envelope  fit  to 
the  dots.  When  the  dots  in  the  upper  left  hand  corner  (straight 
line  overlaps)  are  present,  the  A  value  in  (B-21)  (which  is  the 
asymptote  of  equation  (B-21))  is  set  to  the  minimum  upper  left 
hand  dot.  The  envelope  curves  are  then  least  square  fit  to  the 
minimum  dot  at  each  value  of w. 
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The  meaning  of  Figure  B-10  is  depicted  in  Figure  B-12.  For  any 
delay  time  greater  than  A,  the  aircraft  pair  having  a  crosstrack 
closing  speed  and  an  alongtrack  sepration  which  places  it  in  the 
cell  in  Figure  B-8  is  assumed  to  overlap  while  still  flying 
straight.  For  delay  times  less  than  A,  the  turn  rate  will 
determine  whether  the  pair  will  overlap.  Therefore,  the  shaded 
region  in  Figure  B-12  represents  those  values  of  t<j  and  w  for 
which  aircraft  pairs  will  overlap  given  that  the  aircraft  pair 
was  initially  on  the  conflict  region  boundary  within  the 
specified  ranges  of  x  and  y. 

It  should  be  mentioned  at  this  point  that  for  a  given  cell  in  x 
and  y,  the  horizontal  overlap  region  in  t<j  and  u>  does  not 
always  have  the  form  shown  in  Figure  B-12.  Although  Figure  B-12 
is  the  most  common  form  there  are  three  other  forms  as  shown  in 
Figure  B-13.  Figure  B-13a  shows  the  case  where  the  minimum  time 
delay  is  for  the  straight  line  collision  only.  Thus,  regardless 
of  the  value  of  oj  ,  the  aircraft  pair  will  overlap  for  time 
delays  longer  than  A.  That  is,  for  any  turn,  as  long  as  the 
delay  is  less  than  A,  the  horizontal  overlap  will  be  avoided. 

For  a  given  range  of  x  and  y  values  it  is  easy  to  check  the 
existence  of  a  straight  line  overlap.  It  may  turn  out  that 
there  are  no  straight  line  overlaps  possible  within  the  given 
range  of  x  and  y.  In  this  situation  there  can  be  two  forms  of 
the  tj,  horizontal  overlap  region  as  shown  in  Figures  B-13b 
and  B-13c.  In  Figure  B-13b  the  envelope  of  the  overlap  region 
has  the  same  functional  form  as  the  region  in  Figure  B-12, 
namely  that  of  equation  B-21.  But  since  no  straight  line 
collisions  are  possible,  equation  (B-21)  is  fit  finding  the 
value  of  A  via  least  squares  fitting  of  (B-21)  to  the  maximum 
and  minimum  points  (at  each  w )  in  the  region. 

The  second  case  where  no  straight  line  overlaps  are  possible  is 
shown  in  Figure  B-13c.  In  this  case  there  is  apparently  no 
asymptotic  behavior  to  the  horizontal  overlap  region.  Hence, 
equation  (B-21)  does  not  hold.  In  this  case  the  envelope  of  the 
overlap  region  was  fit  to  the  functional  form 


tj  *  Cln  ui  +  A. 


( B— 2  2 ) 


The  lower  boundary  of  the  horizontal  overlap  region  was  found  by 
least  squares  fitting  the  minimum  values  of  t^  for  each  w. 
Using  the  same  value  for  C,  the  upper  boundary  was  found  by  a 
least  square  fit  on  the  maximum  values  at  each  .  The  same 
value  of  C  was  used  for  the  upper  and  lower  boundaries  so  that 
the  upper  and  lower  boundaries  would  not  cross. 
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B5.  COMPUTATION  OF  THE  PROBABILITY  OF  HORIZONTAL  OVERLAP 

We  now  have  all  the  elements  necessary  to  compute  the 
probability  of  horizontal  overlap.  If  an  aircraft  pair  is 
initially  on  the  conflict  region  boundary  and  has  a  specific- 
time  delay  followed  by  a  specific  turn  rate  then  there  will  be 
overlap  in  the  horizontal  plane.  If  being  on  the  conflict 
region  boundary  is  independent  of  the  time  delay  and  the  choice 
of  the  turn  rate,  and  the  choice  of  the  turn  rate  is  independent 
of  the  time  delay,  then  the  probability  of  horizontal  overlap 
basically  will  be  the  probability  of  being  on  the  conflict 
region  boundary  times  the  probability  of  having  a  particular  set 
of  delay  times  and  turn  rates. 


The  conflict  region  boundary  discussed  in  Appendix  A  is  a 
surface  in  the  three  dimensional  space  x,  y,  and  y  as  shown  in 
Figure  B-14.  To  estimate  the  probability  of  being  on  this 
surface,  the  surface  is  enclosed  with  I  cells  of  width  Ax,  Ay 
and  Ay{,  centered  at  the  points  xi(  yi,  and  yj  for 
i^l, The  probability  of  being  on  the  conflict  boundary  is 
given  by 
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yi\<^yi/2,  |  y  -  y,.|  <^y/2 


As  discussed  previously  in  this  appendix,  the  horizontal  overlap 
regions  (associated  with  the  conflict  boundary)  were  defined  in 
terms  of  cells  in  the  x,  y  plane,  centered  at  the  points  x^, 
and  y{ ,  and  of  fixed  width  Ax  and  Ay.  Given  these  fixed 
values  for  each  i,  there  is  a  unique  value  of  y  (y{),  and  a 
unique  value  of  Ay  ( Ay^)  which  will  enclose  the  conflict 
surface.  These  values  are  derived  from  the  definition  of  the 
conflict  region  given  in  equation  (A-5). 


We  now  turn  over  attention  to  evaluating 
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xi(<Ax/2 


-  yii<Ayi/2 ,  ly  -  yti<^y/: 
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for  any  i.  Now,  the  aircraft  enter  their  respective  routes 
independently.  Therefore,  an  aircraft  pair's  alongtrack 
separation  (before  a  controller  intervention)  will  be 
independent  of  the  pair's  crosatrack  separation  and  crosstrack 
closing  speed.  Thus  (B-24)  ran  he  rewritten  as 

P[CBJ  =  p[|x  '  xi:<4x/2J  *  p[/y  -  yii<  4y./2,|y  -  y/2] 

(  H  -  2  j ) 

The  probability  of  x  (the  alongtrack  separation)  being  in  the 
range  Xj-Ax/2  to  x{  +  Ax/2  will  depend  on  the  traffic  loading 
on  each  route.  The  estimation  of  this  probability  is  discussed 
in  Appendix  E.  As  shown  in  that  appendix  this  probability  is 
constant  over  the  ranges  of  x  under  consideration  (i.e., 
aircraft  pairs  separated  alongtrack  by  less  than  5  nrai).  The 
probability  is  denoted  as  P^x  and  will  depend  on  the  traffic 
loading  and  the  cell  size  Ax,  but  not  the  position  of  the  cell, 
Xj  ,  Thus  (B-25)  can  again  be  rewritten  as 

PCCBi]  =  pAx  *  Ay'  (B-26) 

where 

p(yi»  Ayi,  yi,  Ayj  =  P  |y  -  y  l<Ay  .  / 2,  |  y  -  yj<  Ay/2  j 

The  probablity  P(y[,  Ayi  yi,  Ay)  in  (B-26)  is  the  joint 
probability  that  the  aircraft  pair  has  a  crosstrack  separation 
in  the  ith  cell  and  has  a  crosstrack  closing  speed  in  ith  cell. 
This  joint  probability  can  be  estimated  from  the  data  collected 
by  the  PAA.(^)  Since  the  data  was  gathered  on  individual 
aircraft,  a  convolution  of  the  data  is  required  to  generate  a 
bivariate  histogram  of  crosstrack  separation  and  crosstrack 
closing  speed.  Appendix  F  discusses  the  procedure  used  to 
construct  the  bivariate  histogram. 

For  each  cell  as  shown  in  Figure  B-14  there  will  be  a  horizontal 
overlap  region  in  tj-^j  space.  Such  a  horizontal  overlap 
region  is  shown  in  Figure  B— 15.  From  sections  B2 . 3  and  B3.0  of 
this  appendix  we  have  a  histogram  of  delay  time  and  a  histogram 
of  turn  rate.  These  histograms  are  shown  along  their  respective 
axes  in  Figure  B-15.  Assume  that  there  are  m  cells  in  the  delay 


DELAY 


FIGURE  B-15 

HORIZONTAL  OVERLAP  REGION 
FOR  CONFLICT  BOUNDARY  CELL  i 


histogram  and  n  cells  in  the  turn  rate  histogram.  Consider  the 
kth  cell  in  the  delay  histogram  and  the  t'th  celt  in  the  turn 
rate  histogram.  As  shown  in  Figure  B- 1 1  thi9  combination  of  the 

kth  delay  and  the  fth  turn  rate  is  in  the  horizontal  overlap 

region.  Let  us  define  the  quantity  rj<t'£  with  the  following 
property: 

0  <  rk^iil  and  is  in  the  same  proportion  that  the  ki  cell 
i9  enclosed  within  the  horizontal  overlap  region 

corresponding  to  the  ith  conflict  regon  boundary  cell. 

Let  us  also  define  the  probability  that  time  delay  has  a  value 
in  the  kth  cell  as  P^i  and  the  probability  that  the  turn  rate 
has  a  value  in  the  £th  cell  as  P t.  It  should  be  noted  that 
both  r^^  and  Pk£  depend  on  i,  the  conflict  region  boundary 
cell.  The  values  of  rj^j  will  obviously  depend  on  the 

location  of  the  aircraft  pair  on  the  conflict  region  boundary 
because  there  is  a  different  horizontal  overlap  region  for  each 
position  on  the  conflict  region  boundary.  The  delay  time 

histogram  will  also  depend  on  the  position  on  the  conflict 
region  boundary  because  the  delay  due  to  detection  of  the 

potential  conflict  will  depend  on  the  crosstrack  closing  speed 
of  the  aircraft  pair. 

We  can  now  write  down  the  probability  that  an  aircraft  pair  will 
be  in  horizontal  overlap.  That  probability  is 

I  n  m 

PH  ■  pAx  X  p(v  "V  h-  iy)  X  X  r*kW/  (“-27' 

i-1  Z-l  k=l 

All  of  the  terms  in  this  expression  are  computable  based  on  the 
analyses  in  this  appendix  and  Appendices  C,D,E,  and  F. 
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APPENDIX  C 

PROBABILITY  OF  OBSERVING  AN  AIRCRAFT 
PAIR  WITHIN  THE  CONFLICT  REGION 


Cl.  INTRODUCTION 


This  appendix  addresses  the  applicable  conflict  region,  the 
tracker  errors,  and  the  computation  of  the  probability  of  being 
detected  in  the  confli  t  region. 

Consider  an  aircraft  pair  separated  by  a  crosstrack  distance  y 
and  closing  with  a  crosstrack  speed  y.  What  is  the  probability 
that  the  aircraft  pair  is  detected  to  be  in  potential  conflict 
based  on  the  information  from  one  radar  scan?  The  probability  of 
observing  an  aircraft  pair  as  being  in  potential  conflict 
depends  on  two  factors — the  radar/ tracker  performance  and  the 
conflict  region.  Tracker  performance  is  based  on  the 
propagation  of  the  surveillance  errors  through  the  tracker.  In 
our  model  it  is  assumed  that  the  the  crosstrack  closing  speed 
and  the  crosstrack  separation  error  distribution  is  a  bivariate 
normal  with  zero  mean. 

C2 .  THE  CONFLICT  REGION 


As  defined  in  Appendix  A,  the  conflict  region  is  that  region  in 
the  crosstrack  closing  speed  -  crosstrack  separation  space  where 
an  aircraft  pair,  if  projected  ahead  along  a  straight  path  for  2 
minutes,  will  come  within  5  nmi  of  each  other.  The  shape  of 
that  region  for  the  true  crosstrack  separation  and  crosstrack 
closing  speed  is  shown  as  the  solid  line  in  Figure  C-l. 

In  Section  C4,  we  will  discuss  estimating  the  probability  that 
an  aircraft  pair  is  observed  inside  the  conflict  region  given 
that  the  pair  is  truly  inside  the  conflict  region.  In  that 
estimation  process  we  will  integrate  a  bivariate  normal 
distribution  over  a  polygon.  Thus,  we  have  limited  the  extent 
of  the  conflict  region  by  considering  the  maximum  crosstrack 
closing  speed  to  be  10vymax'  the  maximum  crosstrack 

closing  speed  based  on  data  is  2Vymax,  then  the  choice  of 
lOVymax  should  be  ample  to  include  the  errors  introduced  by 
radar  and  tracker  which  would  indicate  that  the  aircraft  pair  is 
in  potential  conflict. 
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C3.  THE  TRACKER  PERFORMANCE 


The  performance  of  the  tracker  is  characterized  by  errors  in  the 
crosstrack  separation  and  crosstrack  closing  speed.  Those 
errors  are  assumed  to  he  correlated  and  normally  distributed. 
The  manner  in  which  the  parameters  of  the  bivariate  normal  were 
estimated  was  through  a  simulation  of  the  NAS  tracker. 

C3.1  The  Tracker  Simulation 


The  simulation  was  constructed  to  represent  the  performance  of 
the  NAS  tracker  for  discrete  targets  from  a  single  radar.  The 
exact  scenario  which  was  run  was  idealized  for  ease  of  traffic 
generation  but  at  the  same  time  would  show  the  performance  of 
the  tracker  to  be  worse  than  one  would  expect  from  a  typical 
pair  of  aircraft  flying  the  routes. 

The  simulation  consisted  of  flying  a  pair  of  aircraft  on  a  set 
of  parallel  routes.  The  radar  which  observed  these  aircraft  was 
situated  between  and  at  one  end  of  the  routes.  Since  the  radar 
errors  in  azimuth  are  greater  than  the  radar  errors  in  distance, 
this  orientation  of  the  radar  to  the  routes  is  the  worst  for 
estimating  the  crosstrack  positions  and  crosstrack  velocities  of 
the  aircraft.  This  is  another  conservative  aspect  of  this  model. 

Since  the  spatial  relationship  of  the  aircraft  to  the  radar  is 
important,  the  aircraft  were  placed  on  their  respective  routes 
in  alongtrack  proximity  to  each  other.  Alongtrack  proximity  is 
defined  as  being  within  an  alongtrack  distance  for  which  a 
potential  conflict  is  possible.  (In  the  case  of  the  conflict 
region  being  used  this  means  that  the  aircraft  has  to  be 

separated  alongtrack  (interroute)  by  less  than  5  nmi.)  Thus, 
once  the  first  aircraft  was  placed  on  its  route,  the  second 

aircraft  was  placed  randomly  on  its  route  within  alongtrack 
proximity  of  the  first  aircraft. 

The  tracks  of  the  aircraft  were  idealized  to  sinusoids.  When  an 
aircraft  flies  with  turns  (such  as  a  sinusoidal  pattern)  the 
tracker  will  lag  the  aircraft  through  the  turns.  Therefore,  it 
is  important  that  the  sinusoids  of  the  two  aircraft  not  always 

be  in  the  same  synchronization  or  the  results  will  be  biased. 

The  basic  assumption  is  that  the  navigation  on  the  two  routes  is 
independent.  A  particular  synchronizat ion  of  the  sinusoids  for 
each  replication  of  the  simulation  would  definitely  invalidate 
the  assmption  of  independence.  Therefore,  for  each  replication 
the  initial  heading  of  each  aircraft  was  randomly  chosen  from 
among  those  headings  possible  for  the  specified  sinusoid. 
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C3.2  The  NAS  Tracker 


The  NAS  tracker  L9  a  bimodal  alpha-beta  tracker.  The  equations 
which  describe  the  tracker  are  as  follows: 

xsn  3  X  p(n-l)  +  °  (Xrn-Xp(n_i ) ) 

^an  3  ^p(n-l)  *  °  ^rn“^p(n-l 1 ^ 

^xn  3  ^x(n-l)  *  P  iXrn-Xp(n-l ) ) 

^yn  3  ^y(n-l)  ♦  P  ^rn“Tp(n_i))/t 

(C-l) 

xpn  3  ^sn  +  ^xnfc 

^rn  3  Ysn  +  Vynt 


where 


*sn,  ^sn  3  X,Y  position  estimates  for  scan  n 

*pn,  ^pn  3  X,Y  predicted  postion  estimates  made  at 
scan  n  Tor  the  position  at  scan  n+1 

Xrn»  Yrn  3  reported  positions  from  the 

surveillance  system  at  scan  n 


Vxn,  vyn  3  *,Y  velocity  estimates  for  scan  n 
°  *  tracker  position  gain 
P  =  tracker  velocity  gain 

The  bimodal  aspect  of  the  tracker  comes  into  play  through  the 
choice  of  the  parameters  a  and  P.  If  the  predicted  position 
from  the  previous  scan  is  within  a  circle  of  a  given  radius  of 
the  reported  position,  then  it  is  assumed  that  the  aircraft  is 
flying  straight  and  the  appropriate  a  and  P  are  used.  If, 
however,  the  predicted  position  is  outside  the  circle  but  within 
a  circle  of  larger  radius,  then  there  is  a  possibility  that  the 
aircraft  is  making  a  turn  and  another  o  and  P  are  chosen.  If 
the  predicted  position  is  outside  the  larger  circle,  the  return 
is  no  longer  associated  with  the  track  and  the  track  goes  into 
the  coast  mode.  In  the  coast  mode  the  track  is  extrapolated  to 
the  next  radar  update  time  using  the  last  predicted  position  and 
the  last  estimated  velocity.  After  the  track  coasts  for  several 
consecutive  scans,  the  track  is  dropped.  The  values  for  o  and  P 
are  given  in  Table  C-l. 


NAS  TRACKER  PARAMETERS 


SHAFT  SEAFCF  A  DFA  I  TARGE  SEAFCH  AP^A 

(1  NMT  P  ADT  OS)  I  (6  "WI  RADIOS) 


C3.3  Simulation  Results 


The  statistics  of  interest  .-ire  the  mean,  standard  deviation  and 
correlation  coefficient  of  the  crosstrack  separation  and 
crosstrack  closing  speed  errors.  Since  we  are  interested  in  the 
separation/closing  speed  errors  on  a  per  update  basis,  the 
averages  and  standard  deviations  are  taken  over  all  the 
individual  updates  for  all  replications.  This  procedure  will 
give  larger  standard  deviations  than  if,  for  instance,  we  had 
taken  the  mean  separation  and  closing  speed  errors  of  each 
replication  and  then  found  the  standard  deviation  of  these 
means.  The  manner  in  which  the  means  and  standard  deviations 
were  computed  is,  therefore,  more  appropriate  to  the  use  of  the 
stat ist ic  s. 

The  beacon  radar  was  assumed  to  have  a  range  quantization  of 
0.125  nmi  and  an  avionics  bias  error  of  .7  (Jisec  (.113  nmi).  The 
one  sigma  azimuthal  error  of  the  beacon  was  assumed  to  be  .26 
degrees.  If  these  errors  are  assumed  to  have  a  fixed  value, 
then  the  important  parameters  in  this  simulation  are  the  period 
and  amplitude  of  the  sinusoid  pattern  that  the  aircraft  fly. 
Table  C-2  shows  the  variation  in  the  crosstrack  separation  and 
crosstrack  closing  speed  errors  for  three  different 
period/amplitude  combinations.  The  data  that  was  taken  in  the 
Cleveland  ARTCC  has  shown  that  it  is  not  unusual  for  aircraft  to 
wander  4  nmi  off  their  assigned  route  centerline.  The  tracks 
also  oscillated  about  the  route  centerline  and,  while  the  tracks 
were  not  truly  sinusoidal,  they  appear  to  have  periods 
approaching  100  nmi.  A  reasonable  range  of  values  for  the 
tracker  performance  would  be  a  crosstrack  separation  error 
standard  deviation  between  .45  and  .70  nmi,  a  crosstrack  closing 
speed  error  standard  deviation  between  50  and  160  kts,  and  a 
correlation  coefficient  of  -0.8.  For  the  risk  analysis  reported 
in  volume  I  of  this  report  a  value  of  .7  nmi  was  used  for  the 
one  sigma  value  of  the  crosstrack  separation  error,  160  kts  was 
used  for  the  one  sigma  value  of  the  crosstrack  closing  speed 

error,  and  -.8  was  used  for  the  correlation  coefficient. 

The  histograms  of  the  marginal  distribution  of  crosstrack 
separation  and  crosstrack  closing  speed  errors  are  shown  in 
Figures  C-2  and  C-3,  respectively  for  the  case  of  a  period  of  60 
nmi  and  an  amplitude  of  4  nmi,  also  shown  in  these  figures  are 
the  normal  distribution  fits  to  the  simulation  results.  As  one 
can  see  from  Figure  C-2  the  normal  assumption  is  quite  good  for 
the  crosstrack  separation  error  distribution.  The  crosstrack 
closing  speed  errors  are  unimodal  and  symmetrically  distributed 

(see  Figure  C-3).  However,  this  error  distribution  is  not  fit 
very  well  by  a  normal  curve.  The  bulges  in  the  histogram 

between  150  and  300  kts  are  probably  due  to  the  fact  that  the 

aircraft  are  flying  a  sinusoidal  pattern. 
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TABLE 


FIGURE  C-2 

MARGINAL  DISTRIBUTION  OF  POSITION  ERROR 


C4.  INTEGRATION  OF  A  BIVARIATE  NORMAL  DISTRIBUTION  OVER  A  PO LYG ON 

The  problem  of  computing  the  probability  of  observing  an 
aircraft  pair  within  the  conflict  has  been  reduced  to  finding 
the  volume  under  a  bivariate  normal  distribution  over  a 
polygon.  To  illustrate  this,  consider  Figure  C-4.  The  polygon 
defined  by  the  points  A',B',C',  and  D'  is  the  conflict  region 
discussed  above.  Let  an  aircraft  pair  have  a  crosstrack 
separation  of  py  and  a  crosstrack  closing  speed  of  p£  as 
shown  in  Figure  C-4.  The  output  of  the  tracker  will  have  errors 
which  are  assumed  to  be  distributed  normally  in  crosstrack 
separation  and  crosstrack  closing  speed.  The  ellipses  shown  in 
Figure  C-4  represent  constant  probability  lines  from  the 
bivariate  normal  distribution  of  errors  of  the  tracker. 


If  the  polygon  in  Figure  C-4  is  denoted  as  B,  the  probability  of 
observing  the  aircraft  pair  in  that  polygon  given  that  the  pair 
is  at  py  and  p£  is 
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where  is  the  standard  deviation  of  the  tracker's  crosstrack 

separation  error  and  a y  is  the  standard  deviation  of  the 
tracker's  crosstrack  closing  speed  error.  The  quantity  P  is  the 
correlation  coefficient  between  the  errors  in  the  crosstrack 
closing  speed  and  the  crosstrack  separation. 

The  first  step  in  computing  the  integral  (C-2)  is  to  transform 
the  y  and  y  axes  using  the  following  transformation: 
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A 


CROSSTRACK  CLOSING  SPEED 


y-\  +  y 


(C-3) 


v  =  -1 


where 


Making  this  transformation  circularizes  the  bivariate  normal 
probability  density  function  under  the  integral  in  (C-2l.  This 
makes  the  correlation  between  the  random  variables  u  and  v  equal 
to  zero.  Another  important  feature  of  this  transformation  is 
that  straight  lines  are  mapped  into  straight  lines  and  convex 
polygons  are  mapped  into  convex  polygons.  Schematically,  Figure 
C-5  shows  the  result  of  mapping  the  polygon  B  in  y,y  space  into 
u,v  space. 

Since  a  polygon  can  always  be  defined  by  a  set  of  triangles,  we 
will  compute  the  probability  over  the  polygon  by  accounting  for 
the  probability  over  a  set  of  triangles  that  define  the 
polygon.  Therefore,  consider  a  triangular  area  in  u,  v  space  as 
shown  by  the  triangle  OA'B’  in  the  Figure  C~6a/b.  The 
probability  of  being  within  that  triangle  can  be  computed  by 
standard  methods.  Let  the  point  S  be  on  the  line  connecting  A' 
and  B',  such  that  SO  is  perpendicular  to  the  line  A',  B' .  If 
the  point  S  lies  between  B'  and  A',  as  shown  in  Figure  C-6a  then 
the  probability  of  being  within  the  shaded  triangle  is 


V(h,k,)+V(h,k2) 


(C-4) 


where 


h  is  the  distance  from  the  origin  to  point  S, 
kj  is  the  distance  from  point  S  to  point  B', 


t 
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THE  TRANSFORMED  TRACKER  OUTPUT  ERRORS  AND  CONFLICT  REGION 


k2  is  the  distance  from  point  S  to  point  A',  and  the 
function  V  will  he  defined  helow. 

If  the  point  S  does  not  lie  between  B'  and  A',  as  shown  in 
Figure  C-6bf  then  the  probability  is  given  by 

V(h,k2)  -  V(h,k!)  (c_5) 

where  k2  >kj 

As  derived  in  Reference  (.1),  the  function  V  can  be  written  in  a 
Taylor  expansion: 


OH)'1  [  >  ( l  *-')  -U’ 


(r-6) 


where  K=k/h<l,  a=l/2h2  and  h>0  and  k  ->  0.  If\>l,  then 
equation  (C-7)  can  be  used  in  conjunction  with  (C-6): 

V(h,k)  +  V(k,h)  =  (  <P  (h)— 1/2)  ( 0 (k)-l/2 )  (C-7) 


^(h) y  exP(-t2/2)dt  (c_8) 

-<x> 

Returning  to  Figure  C-b,  we  can  see  that  the  values  of  h,  kj  f 
and  k2  can  be  determined  from  the  following  expressions: 


1U2V1  ~  ulv2 I 

^2"^  +  (v2"vjf 


k 


2 


where  b'  3  j  , 


(C-9) 
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To  compute  the  probability  over  the  polygon  shown  in  Figure  C-S, 
we  could  divide  the  polygon  into  triangles  as  shown  in  Figure 
C-7.  In  this  particular  example  the  aircraft  pair  is  within  the 
conflict  boundary  so  that  the  origin  in  u,v  space  will  be 
within  the  transformed  conflict  region  polygon.  This  means  that 
the  probabilities  of  being  within  triangles  I,  II,  III,  and  IV 
in  Figure  C-7  are  added.  If  the  aircraft  pair  were  not  within 
the  conflict  region  then  the  u,v  origin  would  not  be  within  the 
transformed  conflict  region  polygon.  Figure  C-8  shows  such  a 
configuration.  To  compute  the  probability  of  being  within 
polygon  A'B'C'D'  one  would  find  the  probabilities  of  being 
within  triangles  I,  II,  III,  and  IV  and  then  compute  the 
probability  as 

P(A' B'C'D' )=*P( I)+P( II )-P( III)-P(IV).  (C-10) 
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THE  PARTITIONED  TRANSFORMED  CONFLICT 
REGION  (ORIGIN  INCLUDED) 
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APPKNDIXJ) 

THE  FAST  FINITE  FOUR  IKK  TRANSFORM 


1)1.  INTRODUCTION 

The  purpose  of  this  appendix  is:  1)  to  give  a  general  des¬ 
cription  of  some  of  the  properties  of  the  finite  Fourier 
transform,  including  the  properties  which  enable  the  "fast" 
calculation  of  the  transform;  2)  to  describe  the  use  of  a 
fast  Fourier  transform  computer  program  for  computing 
convolutions  and  3)  to  present  a  fast  finite  Fourier  transform 
computer  program.  There  are  many  applications  for  the  fast 
finite  Fourier  transform  which  are  not  treated  here,  such  as  for 
the  estimation  of  spectral  density  functions  and  Fourier  series. 
Articles  in  Reference  9  include  a  discussion  of  the  fast 
Fourier  transform  algorithm  and  its  history,  as  well  as  a 
treatment  of  some  applications. 

The  finite  Fourier  transform,  F(n),  of  a  finite,  complex  series, 
fft)  where  t=0,...,T-l  is  given  by  the  series 

Frnf  =  1  V  2;rint 

"s'?  f(t)  e  T  for  r>=o,  ■  •  -,T-1 .  (D-l) 

The  finite  Fourier  transform  is  the  finite  analogue  of  the  usual 
discrete  Fourier  transform  and  enjoys  many  of  the  properties 
associated  with  the  infinite  dimensional  version.  The  pro¬ 
perties  of  the  finite  Fourier  transform  are  discussed  in  Section 
2  of  this  Appendix.  The  fast  Fourier  transform  is  a  computation¬ 
al  procedure  used  for  efficently  calculating  the  finite  Fourier 
transform.  It  is  not  in  estimation  procedure,  but  simply  an 
efficient  calculation  procedure.  The  fast  Fourier  transform 
algorithm  permits  calculation  of  (D-l)  in  approximately  Tlog^T 
operations  instead  of  the  T^  operations  required  for  its  direct 
calculation.  The  properties  of  (D-l)  which  lead  to  this  effi¬ 
cient  calculation  procedure  are  outlined  in  Section  3  of  this 
appendix.  Section  4  of  this  appendix  discusses  certain  features 
of  the  finite  transform  which  are  necessary  for  the  correct  use 
of  the  fast  Fourier  transform  computer  program  and  for  the  correct 
interpretation  of  the  results  when  it  is  used  to  calculate  con¬ 
volutions.  Finally  one  computer  program  which  can  he  used  for 
calculating  the  fast  Fourier  transform  is  presented  in  Section 
3. 

D2.  PROPERTIES  OF  THE  FINITE  TRANSFORM 

First,  since  the  finite  Fourier  transform  defined  hv  (D-l) 


D-l 


is  a  finite  sum  of  finite  valued  I une t ions ,  it  always  exists 
and  is  finite.  The  proper!  ies  presented  in  this  sort  inn  do  not 
require  any  assumptions  (other  than  finiteness)  roneern  i  up. 
the  function  f(t).  However,  restrictions  on  f(t)  are  fre¬ 
quently  required  for  estimation  procedures  based  on  the  linite 
t  ransf orm. 


If  the  finite  Fourier  transform  of  a  function  f(t)  is  defined 
by  equation  (D-l),  then  us  i  ng  the  fact  that 


T  -  1  2tt  i  nt  -2  n  i mt 


t  =  0 


T  if  n  ^  m  mod  T 
0  otherwi so. 


<  !>-2) 


it  is  possible  to  show  that  the  inverse  relationship 
j.  T-l  -2  Hint 

f(t)  =  -  ^  F(n)  '  T  for  t=0 _ _  T-l  (0-3) 

n=0 

also  holds.  Thus,  the  finite  Fourier  transform  satisfies 
conditions  which  are  similar  to  the  infinite  dimensional 
discrete  Fourier  transform.  Functions  such  as  f(t)  and  F(n) 
which  satisfy  ( I >— 1 )  and  (0-3)  are  referred  to  as  finite  trans¬ 
form  pairs. 


15 v  definition  of  the  finite  transform  in  (D-l),  it  is  ease 
to  see  that  F(n)  is  periodic  with  period  T.  Since  e~,T''  * 
for  any  integer  k,  wo  have  that 


1  T-’ 
F(n+kT)  -  ^ 

I  (t  )  e 

2 it  i  (n+kT)  t 
T 

■  T  t-d 

T-l 

2rr  i  n  t  3 

=  I  V 

1  (t  )  e 

’Y" 

•  T  t=() 

w 

Mi 

2r  \nl 

f  ( t  )  e 

T  * 

\  T  t=0 


for  any  integer  k  and  n=t>, .  .  .  ,  T-l  . 


Bv  the  same  argument,  equation  (I)  ()  implies  that  f(t)  = 

f ( kT+t )  for  any  integer  k  and  t=u . T- I .  Tims,  the  finite 

transforms  are  computed  under  the  assumpt ien  that  fit)  is 
periodic  with  period  T. 

A  useful  property  of  tin1  finite  transform  is  that  if  f(t)  is 
a  real  valued  function,  then 

F(n)  -  f*(T  -  n)  (p-4) 

for  n  =  l i . T/ 2 ,  where  *  denotes  complex  conjugate.  Thus, 

F(n)  is  uniquely  determined  bv  its  values  at  the  points 
n  =  0,  1,  .  .  . , T / 2 . 

Finally,  the  most  important  property  for  the  analysis  pre¬ 
sented  in  this  paper  is  that  convolution  and  point  bv  point 
multiplication  are  dual  operations  in  time  and  frequence,  i.*-., 
if  f(t)  and  F(n)  are  finite  transform  pairs  and  g(t)  and 
G(n)  are  finite  transform  pairs  with  the  same  domain  • , , t  = 

0,  ■  .  .  , T-l ,  tlien 

T  -  1  T-l  -2"hu 

S  f(t  -  u)  r(u)  =  V  F ( n )  C( n)  c  1  .  (p_%) 

U  =  0  0 

Because  of  the  operational  efficiency  of  the  fast  Fourier 
transform  algorithm,  the  calculation  of  the  left  liand  sub 
of  (D-5) proceeds  most  efficiently  by  calculating  F(n)  and 
C(n)  via  (D-I)  and  then  using  (!)-))  to  obtain  the  inverse 
transform  of  F(n)  0(n).  This  method  will  he  most  efficient 
as  long  as  the  total  number  of  operations*  required  to 
calculate  the  left  hand  side  ef  f  D—  S )  direct  Iv  (T-* )  is 
greater  than  the  number  required  to  calculate  F(n),  C.(n)  ami 
the  inverse  transform  of  F(n)  0(n)  OTIogfl).  Thus  if  T  is 
a  power  of  2,  using  the  fast  Fourier  transform  as  described 
above  to  calculate  the  left  hand  side  of  (!>-s)  is  more 
efficient  than  direct  calculation  as  long  as  T?16. 

THE  DOUBLING  ALGORITHM 

Consider  a  series  f(t)  of  length  T  =  2N.  Then  the  finite 
transform  of  f(t)  is  given  by  (1)-1  )  to  be 


* 


Complex  multiplications  and  adds 


V  I 


I  ,  viL'ii.. 

I'  (  ii  )  ■  -  1(1)  ■  2N  for  a  -  <  i ,  •  •  • ,  2N- 1  .  ( |)-h  ) 

'-2N 


hot 


f  (t) 

() 

=  f (2t+l  )  lor 

1=0, 

>  the  series  ol 

odil  numhe n  il  < 

I ement s  ol  f  (t ) 

and 

ft.(n 

:  f (2t  )  lor  t  - 

<),•••,  N  -  1 

he  t  h-  series  of 

even  numbered 

el ement s  of  f (t ) 

Then  the 

finite  transforms 

of  fit)  and 

ft,(i)  are  given 

by  (D-l)  to  be 

1 

N-l 

J  n  i  •  i : 

F0(n)  -  i 
•  N 

21  '»(t)  - 

N 

t=  U 

anti 

Ft,(n)  -  - 
' .  N 

N-l 

int 

E  f-(u  - 

N 

t=n 

Now,  from  ( f)—  f> ) 

we  eat.  write  F(n)  in  terms  of 

F  ,  (n)  and  F  ( 

as  foil ows 

N-l 

2tt  i  u(2t  +  l ) 

N-l  2- 

2N  F(n)  - 

£  1  ( 2 1  1  ) 

e  2N  + 

£  f  (2t)  e  ’ 

t=d 

or 

.  .  N-l 

2 r  int 

2r  i  it  , 

\2F(n)  =  e 

■  2»  A  £ 

f  (t )  e  N 
o 

N 

t  =  11 

N-l 

+  :  E 

2  rr  i  n  t 

f  ( t  )  e 

for  n=l>,  •  •  •  ,  2N  . 

I>-4 


Thus,  since 


2  n i (n+N ) 
2N 


— t* 


i  n 


t  ur  n=<),  •  •  ■  ,  N-  (  , 


e 


2N 


we  have 


*.2F  ( n ) 


2rrin 

2N!  F  (n)  +  F(,Cn)  for  11=0,.  .  .  <  M—  1 


F(|  (n)  4-  Ft,(n)  for  n=N,  •  •  2N-1 


V  (l»-7) 


Tf.e  c-cvat  ion-^  (\}-l )  arc  tin*  so-called  doubling  algorithm.  (avcn 
the  finite  transforms  of  the  two  series  consisting  of  the  odd 
and  even  numbered  elements  of  a  series  of  length  T=2N,  (D- < ) 
•rabies  calculation  of  the  transform  ot  the  total  series  be 

complex  multiplications.  Thus,  if  we  denote  by  the  number 
of  multiplications  required  for  calculating  tlu*  transform  of 
an  .M  point  array,  using  (0-7)  to  calculate  the  transform  of 
a  2M  point  array  wo  have  that 


2.M 


2u(-  +  . 


S  i ns <■  <K  =  1  ,  t  li  i  •  implies  [  ha  t 


r  i>t 


■qN'  log  ,.' 


Thus,  the  fast  Fourier  transform 
of  complex  multiplications  from 
calculation  to  ^Tlog^T.* 


algorithm  reduces 
9 

the  T“  required  bv 


the  number 
the  direct 


If  the  length  of  the  data,  T,  is  a  power  of  two,  the  fast 
Fourier  transform  proceeds  as  follows: 


1,  The  series  of  data  is  split  into  I  subseries  of 
length  1.  Thus,  the  Fourier  transform  for  each  series 
is  the  series  itself.  These  T  transforms  are  called 
the  first  level  transforms, 

1  The  figure  Tlog?T  given  on  page  D- 1  is  obtained  by  doubling  ^Tlog^T 
in  order  to  include  the  number  of  complex  additions. 
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2.  'Mu*  doub  I  |  ng  .  I  I  ;•(  >  I  i  i  |  mi  III-/)  |  .;  1 1 1 1 .  i  j  ,  I  iii  l  lii  ' 

■  mil  I  /  2+k*  h  I  i  r  s  t  1 1'  vi  I  t  t  i n >.  I  i  ■  t  ins  i  1 1  i  i  i!ii|> 1 1 1 1 ■  [lu  !  i  u  !  i  . 
transforms  for  t  lie  uni  I .  •  v  ■  1  subset-l.-s  k,  where 
k=  0,  •  ■  ■  , T  /  2 ,  Tlicse  arc  i  In  I  i  n  i  1 1-  t  r.ms  t  .  ■  pii . •  t  t  ! u 
17  2  subserics  of  1 1 •  i i > ■. t  ) ■  2, 

1.  Doubling  mitimu.>s  in  this  wav  until  final  lv  t  lu 
transforms  of  the  two  subscri  rs  of  length  T/2  art 
combined  via  tho  doubling  ilgurithw  to  give  the 
finite  trail  storm  of  the  enliri  series. 

The  fast  Fourier  transform  program  presented  in  Section  ’>  opor.Ui 
solely  on  data  arrays  of  length  a  power  of  two.  Other  com¬ 
puter  programs  (using  similar  algorithms)  can  handle  seru'- 
of  any  length.  In  general,  the  restriction  to  a  length  et  a 
power  of  two  has  not  been  Luc  restrictive.  The  series  is 
sir.pl  •  filled  with  zeros  to  give  it  the  proper  length. 

FEATJ  RES 

This  Section  will  discuss  several  features  of  the  use  ot  tin 
fast  Fourier  transform.  These  include  features  relating  to 
the  input  of  data  for  'ransformat  ion,  the  required  length 
of  t  fie  input,  and  the  interpretation  of  the  output  when  tin- 
fast  Fourier  transform  is  used  to  perform  convolutions.  Knowl¬ 
edge  of  these  features  was  ne. essarv  for  correct  use  of  tin- 
fast  Fourier  transform  in  the  anal  vs i s  presented  in  this 
paper.  A  more  complete  description  of  the  fast  Fourier 
transform  and  features  applicable  in  other  types  of  analysis 
mav  be  found  in  Heferen<i-9. 

As  mentioned  previouslv,  the  finite  Fourier  transform  one rates 
under  the  assumption  that  t he  function  to  In  transformed  is 
periodic  with  period  T.  Thus,  the  a  Igor  it  lint  calculates  tin- 
finite  transform  of  a  periodic  extension  of  the  data  observed 
between  t=h  and  L=T-I.  If  the  data  :(t)  represents  observa¬ 
tions  which  are  symmetric  about  zero,  i.e. ,  if  actual  observa¬ 
tions  were  taken  at  t.  =  -  1  /.:+  I  ,  •  •  • ,  •'/ j  then  the  data  should 
be  shifted  in  the  following  wav  before  its  finite  transform 
is  calculated.  Define 

:<  ( t  )  for  t r  1 J ,  .  .  •  ,  1/- 

y  (t  )  =  •  (  D  -  8) 

x ( t - T )  lor  t  =  1  /  2  +  1,  .  .  •  »  I-l. 


l)-t 


Then  the  .U-s  t  red  t  r,iiisf»rm  of  x(t)  i  .  obi  a  ined  lc:  trans¬ 
forming  •:  ( t  ).  ti  would  also  ho  possible  to  transform  I  ho 
function  vfu  -  x ( t  -  J/,  +1)  tor  t  =  0,  •  •  •,  T- 1  instead.. 

This  represents  a  shifted  version  of  tlio  poriodio  oxtonsion 
c>f  (D-b )  and  its  transform  would  just  ho  a  shifted  version 
of  the  transform  of  (|)-8f.  (D-8),  however,  directly  repre- 

sonts  the  function  which  is  symmetric  about  zero. 

There  ire  several  features  of  the  finite  Fourier  transform 
which  must  he  c  onsw'nred  when  it  is  used  to  calculate 
■onvolut  ions.  First  ,  lot  us  ei'nsider  the  convolution  of  two 
funi  t  ion  lit  )  and  p(t  ).  Lot  us  assume  that  f(t)  was  ob¬ 
served  for  t  -  0,  •  •  •  ,  1-1  and  it  ( t  )  was  obs orved  for  t  =  0  ,  •  •  •  ,  N-  1  . 
Let  T  he  a  power  of  two,  with  the  property  that 

i  '  \.  +  N  -  l  (If-9) 

Th--r,  define  f(t)  -  <  for  t  =  ].,•••,  T-i  and  g(t)  =  0  for 
t  -  >!,•••,  T- 1 .  TIk  convolution  of  f(t)  and  g(t)  can  then 
he  computed  via  f D— " ) .  One  way  to  justify  the  constraint 

•  D—  9 >  is  to  view  the  convolution  in  (D-5)  as  the  probability 
function  of  the  sum  >f  twe'  independent  random  variables  with 
density  functions  f ( t  )  and  gft)  respective  y.  Then  if  f(t) 
is  defined  for  I.  data  points  and  g(t)  is  defined  for  N  data 
points,  the  sum  of  t  lie  two  random  variables  will  have  a  domain 
.  f  definition  ot  M+i.-l  data  points. 

1  he  t  inite  Fourier  t  ranslc  rm  presented  above  is  defined  for  a 
function  l(t)  observed  at  values  of  t  =  0 , • -  - ,  T- 1 ,  In  practice, 
observations  are  frequent  Iv  taken  at  a  fixed  sampling  interval, 
say  A ,  wit  i  l  h  is  not  necissarilv  unitv.  In  .addition,  the  first 
observation  of  real  data  mav  not  correspond  to  an  observation 
at  zero.  For  example,  in  the  analysis  presented  in  this  [taper 
it  was  desired  to  compute  the  convolution  of  two  probability 
functions,  symmetric  about  zero,  observed  in  1/2  mi  intervals. 
There  was  no  zero  observation;  instead,  observations  within 
1/2  mi  of  zero  fell  info  one  of  the  intervals  [-.5,0.)  or 
[(!.,.  5).  The  finite  transform  can  still  lie  used  to  calculate 
the  convolution  of  such  functions.  The  following  analysis 
provides  the  kov  to  into  ■  ,  tat  ion  of  the  results. 

•  •  %  r 

Let  us  assume  that  we  have  two  sets  of  T  data  points  x(s) 
and  y(s)  (including  the  necessary  number  of  zeroes  to  satisfy 
(D-9))  observed  in  time  intervals  of  A  starting  at  time  a. 

Thus  x(s)  md  v(s)  are  defined  for  s  =  At  +  a  where 

,  1-1.  In  order  to  calculate  the  convolution  of  x 
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-'inti  y,  we  set 


lit)  (Al 

♦ 

a) 

and 

g ( 1 )  v  (Al 

+ 

■i) 

f 

or  l  <),•■■  1 

(D-K)) 

The  functions  f(t)  and  g( t )  n 

(* 

now 

i  n 

a  form  to  use 

t  lie  fast 

Fourier  transform  methodology 

L  O 

cal 

cul 

ate  tile  convolution.  UY 

have  from  (D-5),  that  the  output 

o  I 

the 

fast  Fourier 

t  rans  form 

convolution  computation  would 

he 

given 

by  h(t),  where 

’ 

T-  1 

h ( t )  =  £  f(t-u)  g(u) 

t=u. 

,  T-l. 

u=i  I 


Substituting  for  f  and  g  from  (D-10),  wo  have 
T-l 

h(t)  =  D 

u=o 

This  can  be  written  as 


T-l  /  v 

xjA(t-u)  +  aj  v(Au 
u=D  '  / 


+  a). 


T-l 


u=  ( ) 


h(t)  =  £  x  ((At  +  2a)  -  (Au  +  a)J  v  (Au  +  a). 


Since  the  expressions  involving  u  in  the  abov  -  equation  are  the 
same,  this  ran  be  written  as 


h(t)  =  z  (At  +  2a)  for  t  =(i ,  -  ■  -  ,  T-l.  (D- 1  1  > 

In  (0-11),  z  (At  +  2a)  represents  the  desired  convolution  of 
x(s)  and  y(s).  (D-ll)  implies  that  the  series  of  points 
representing  the  convolution  of  two  functions  are  observed 
at  the  same  time  interval  as  tla  original  series,  namelvA. 
However  the  fact  that  the  first  observations  in  the  initial 
series  were  offset  from  zero  bv  "a"  units,  implies  that  the 
first  point  in  the  convolved  series  will  be  offset  from 
zero  by  "2a"  units. 

For  the  example  mentioned  above,  the  observation  intervals 
were  1/2  nmi ,  so  that  A  =  1/2.  The  first  observation  intervals 
were  centered  at  1/4,  so  that  a  =  1/4.  Thus,  the  convolution 
of  x(s)  with  y(s)  has  the  interpretation  (D-ll),  or 

-  z  (t  /2  +  . ’>  )  for  t  =  <> ,  •  •  •  ,  T-I  . 


D-8 


hit  ) 


'I  bus,  t  Ilf  t  i  r  • .  L  i  t  .'in  in  l  In-  array  ri'turiifil  hv  t  lie  program, 
namely  li(0),  represents  the  desired  convolved  density 
tini't  ion,  z  irv.rl  ii.it  ed  it  I/.'  nmi.  Because  of  t  lie  periodicity  ol 
the  finite  Fourier  t  ra  ns  t  ■  ■  r  m ,  we  iiave  that  the  last  item  in  the 
data  arrav  lift;,  name  I  v  h  (  I- 1  )  is  equal  to  h(-l).  Thus,  we  have 
that  i('\  i :  )  ■-  -''Ol,  so  th.it  the  last  item  in  the  array  returned 
bv  the  finite  transform  computer  program  represents  the  desired 
convolverl  density  function  evaluated  at  zero. 


The  same  ciiaracterist  ic  applies  when  (It-a)  is  used  to  determine 
the  convolution  of  three  shifted  f uni  t  ions  as  was  done  in 
the  analysis  discussed  in  this  paper  to  include  distributional 
characteristics  of  delay,  bank  angle, and  detection  (Appendix  Cl. 

In  this  case,  it  can  easilv  be  verified  that 

h  ( t  )  -  .-(At  F  '3a )  for  t  =  0,  •  4  •  ,  T-I.  (11-12) 

Let  us  assume  t  hat ,  A  =  1/d  and  a  -  1/4  as  before.  Thus,  the 
first  cell  of  the  arrav  h(t),  (rorrespond ing  to  t  =  0) , 
actually  represents  t  lit*  convolved  dens  it  v  data  .  centered  at 
.75  nmi.  The  last  item  in  the  array  h(l)  (  orrespondi ng  to 
t  =  -1)  represents  the  convolved  density  dila  centered  at 

. 2  j  nmi . 

D5.  _  THK  COMP l TKR_  PROCRAM 

The  fast  Fourier  transform  program  used  in  this  analysis  is  a 
FORTRAN  IV  program  which  was  obtained  from  Dr.  Robert  S hutnway 
of  Oeorge  Washington  University  and  is  attributed  to  Norman  Brenner 
of  MIT  Lincoln  Lab.  It  is  listed  in  Table  D-l.  The  call 
statement  to  the  subroutine  is: 


CALL  I- FT  f  N’N ,  DATA,  SICN'l) 

whe re  , 

DATA  is  a  complex  double  precision  array  of  length 
2^‘  .  It  contains  the  data  to  be  transformed 
on  tile  call  and  the  transformed  data  on 
ret  urn. 


D-'l 


r.Miu  u- 1 


Tilt  l-A'.T  UNITE  FOURIER  TRANSFORM 

computer  pkoukam 


SUBROUTINE  FFTINN .DAT  A ,S  I  UNI  I 
IMPLICIT  REAI.*H  (A-H.w-Z) 

DIMENS ION  DATA  (1) 

DATA  TWOPI/6. 28 31 85 307 2D0/ 

C 

C  THIS  SUBROUTINE  COMPUTES  THE  FAST  UNITE  FOITRITK  TRANSFORM  OF  A 
C  -  ’ **N-N  POINT  COMPLEX  DOUBLE  PRECISION  ARRAY. 

C  I  NT  UTS:  DATA!.  I  COW  LEX  l.»2**NN  POINT  ARRAY 

C  NN  POWER  OF  ’  CORRESPOND!  NO  TU  THE  LENC.TH  OK 

I.  DATA!  .  ) 

C  SICNl  1.  OK  -I.  ( DOUBLE  PRECISION)  DFPENDINC  ON 

C  WHETHER  TOP  DIRECT  OR  INDIRECT  TRANSFORM  IS 

C  DESIRED 

C  OUTPUTS:  DATA  ( .  J  COMPLEX  DOUBLE  PRECISION  1.  POINT  ARRAY  OF 

TRANSFORMATIONS. 

C  DATA(N)  •  SUM  FROM  T  *  0  TO  T  *  L-l  OF 

C  DATA  (T )  *  EXP(SIONI*2*PI*l*N*T/L) 

C  WHERE  N  »  0,1 . l.-l  AND  I  =■  SQRT1-1.) 


N-2**  I.NN+1 ) 

J=1 

DO  5  I-1,N,2 
IFn-J)I,2,2 

1  TEMP  R”  DATA  ( J) 

TEMP  I -DATA  OD 
DATA(J)*OATA<I) 
DATA(J+1)-DATA(I+1 ) 

DATAII  )-TEMPR 

DAT  A I  T-El  )=*TEMPI 

2  M-N/2 

3  IF(J-M)8,5,A 
>  J-M 

M»M/  2 

IFIM-2) 3, J,  ( 

V  I-  J+M 

MMAX*2 

b  I F IMMAX-N) 7,10,10 

7  ISTEP-2*MMAX 

TH  ETA*  S I  GUI  *TW0P  I  /  n.O  AT  (  MM  AX ) 
S I NTH" DS IN (THETA/ 2 . PO) 

WSTP  R—  2  .  D0*S  I  NTH*S  I  NTH 
WSTPI»DS  IN  (THETA) 

VK*1  .DC 
WI-0.00 

DO  9  M-1.MMAX.2 
DO  8  I  "M ,  N ,  I  STEP 
J-t-WKAX 

TEMPR-WR*DArA(J)-WI*DATA( J+l ) 
TEMPI"WR*DATA(J+1)+DATA( J)*WI 
DATA(J)-DATA(I)-TEMPR 
DATAIJ+1  )- DATA  (I -El) -TEMP  I 

data(i)*data(I)+tempr 
2  data(i+1)-data(i+d+tfmpi 

TEMPR"WR 

VR*WR*WSTPR-WT  AWSTPI+WR 
9  VI-W!  *WSTPR*TFMPR*WSTPI+WI 

MMA.X-  ISTEP 
CO  TO  b 
RETURN 
END 


D-l  0 


NN  is  the  power  of  two  corresponding  to  tli 

length  of  the  data  series. 

S1GNI  is  double  precision  +1  depending  on 
whether  a  direct  or  inverse  transform 
is  desired. 
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INTRODUCTION 


PROBABILITY  OF  AI.ONCTRACK  SEPARATION 


El. 


In  Appendix  B  there  is  a  term  P^x  which  is  the  probability  of 
that  two  aircraft,  cleared  on  adjacent  parallel  routes,  have  an 
alongtrack  separation  of  Ax/2.  This  appendix  will  show 
that  P is  constant  over  the  range  of  alongtrack  separations 
being  considered.  This  will  be  followed  by  a  discussion  of  the 
relationship  between  P  ^x  and  Px,  the  probability  of 
alongtrack  overlap.  The  estimation  of  Px  using  an  analytical 
model  will  then  be  addressed. 

E2.  DISTRIBUTION  OF  ALONGTRACK  SEPARATION 

To  show  that  P  is  constant  for  |  x+  Ax/2|<  D  we  first 

investigate  the  distribution  of  alongtrack  separation  between 
aircraft  that  are  assigned  to  adjacent  parallel  routes. 
Consider  the  situation  depicted  in  Figure  E-l.  The  aircraft  on 
the  same  route  are  separated  by  a  distance  W.  The  air  traffic 
control  rules  state  that  W  has  to  be  greater  than  some  minimum 
distance  D  (the  ralar  separation).  Since  W  is  a  random 
variable  it  has  a  piobability  density  furction  which  can  be 
written  as 


f(W) 


o  0< W<D 

g(W)  D<W 


(E-l) 


where  g(W)  is  an  arbitrary  probability  density  function. 

Now,  the  aircraft  on  route  2  are  assumed  to  operate 

independently  from  the  aircraft  on  route  1.  Thus,  if  we  have 
aircraft  A|  and  on  route  1,  then  an  aircraft  A2  on 

route  2  would  in  general  be  positioned  between  the  two  aircraft 
on  route  1  as  shown  in  Figure  E-l.  The  alongtrack  distance  from 
A2  to  A]  is  denoted  as  x-j  while  the  alongtrack  distance 
from  A2  to  B-[  is  denoted  as  X2.  Since  the  position  of 

A2  with  respect  to  A]  andB^  is  randomly  chosen,  the 
distances  x-[  and  X2  are  random  variables. 
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INTERAIRCRAFT  SPACING  GEOMETRY 
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From  renewal  theory  (Reference  10)  we  have  two 
first  result  is  that  the  random  variables  x-j 
identically  distributed.  The  second  result  i 
probability  density  function  is  given  a9 


resul ts . 

and  *2 
s  that 


The 
are 
th  i  s 


h  ( x ) 


(  E-2  1 


But,  we  are  only  interested  in  the  probability  density  function 
of  x  (i.e.,  Xj  and  X2  in  Figure  E-L)  in  the  range  of  0  to 
D.  One  recalls  that  this  is  our  definition  of  proximity  from 
Appendix  A.  Thus,  if  we  are  considering  0  x  D,  then  equation 
(E-2)  becomes 


’n(x) 


x<D 

1  -J  0  dW 

"E(W)  =  ECW) 


(E-3) 


where  E(W)  is  the  expected  value  of  the  variable  W.  Equation 
(E-3)  tells  us  that  the  probability  density  function  of  x  is 
constant  over  the  range  Q<  x<  D.  This  means  that  it  is  equally 
likely  to  find  adjacent  route  aircraft  separated  betwen  0  and  D 
nmi . 

E3 ,  RELATIONSHIP  BETWEEN  PAx  AND  P, 

In  the  notation  of  Appendix  B 

P4X  =  P£|X  "  Xi|<^/2]* 

f  E-4) 

For  Jx£+^x/2|<D,  the  results  of  renewal  theory  (equat'on  (E-3) 
tell  us  that 


P^x  =  Ax/E(W)  . 


(E-5) 


From  the  procedural  collision  risk  model  we  have 


E4 . 


Px  =  P[x<  2  Ax] 


Since  2  X  x  <  D,  and  the  distribution  of  alongtrack  separation 
between  aircraft  on  adjacent  routes  is  uniform  within  the  range 
0  to  D,  we  can  write 


Ax 


(E-7) 


Rearranging  terras  one  arrives  at 


(  E-8 ) 


for  J  +  Ax /2  |  D. 

Therefore,  we  have  determined  that  P^  is  a  constant  that  is 
related  to  Px  as  shown  in  equation  (E-8).  The  next  section 
will  discuss  how  one  can  estimate  Px. 

ESTIMATING  Px 

An  analytical  model  to  estimate  Px  was  developed  by  the  ICAO 
Review  of  the  General  Concept  of  Separation  Panel  (RGCSP) 
(Reference  11).  For  aircraft  that  are  coaltitude,  in  level 
flight,  and  flying  the  same  direction,  the  analytical  model 
gives  a  value  of  Px  as 


Tl+!,2) 


(E-9) 


where  Xx  is  the  length  of  the  aircraft, 

V  is  the  average  velocity  of  the  aircraft,  the  average 
velocity  being  the  same  on  both  routes,  and 


K-4 


is  the  flow  rate  on  route  i. 

The  validity  of  this  analytical  model  was  tested  by  comparing 
the  results  from  (E-9)  with  data  on  observed  aircraft  passings 
in  an  operational  environment.  This  comparison  wa9  made  using 
data  from  a  two  week  period  in  the  FAA'3  large  scale  data 

collection  in  the  Cleveland  ARTCC  (Reference  12). 

Since  equation  (E-9)  was  developed  under  the  assumption  that  the 
flow  rates  are  constant  during  a  steady-state  period,  it 

was  necessary  in  Reference  12  to  determine  the  flow  rates  over  a 
set  of  assumed  steady-state  periods.  Time  periods  for  1/2  hour 
to  6  1/2  hours  were  investigated.  The  ratio  of  the  total  flying 
time  in  passing  to  the  total  flying  time  was  computed  for  each 
day.  These  values  were  compared  to  the  Px  values  calculated 

from  the  model  under  the  various  steady  state  assumptions  for 
each  day  and  for  the  entire  time  period. 

The  results  from  the  investigation  of  the  data  showed  that  the 
interarrival  spacings  are  consistent  with  a  Poisson  process. 
Also,  there  was  a  diurnal  pattern  observed  in  the  data.  From 
day  to  day  there  was  some  variation  in  the  arrival  rates  for  a 
particular  "steady-state"  time  period  but  not  more  than  one 
would  expect  from  a  Poisson  process.  This  result  allowed  the 
combination  of  the  data  in  the  same  "steady-state"  period  across 
all  days.  It  was  concluded  in  Reference  16  that  the  1/2  hour 
"steady-state"  time  period  gave  the  best  approximation  to 
steady-state.  For  this  steady-state  time  period,  the  data 

derived  values  of  Px  and  the  analytical  model  derived  values 
of  Px  were  comparable  over  a  long  time  period  for  same 

direction  flight  on  both  routes.  Therefore,  for  properly 
interpreted  data,  the  analytical  model  provides  a  vehicle  for 
estimating  the  fraction  of  the  flying  time  spent  in  passing. 


APPENDIX  F 

CALCULATION  OF  THE  JOINT  CLOS ! NC  SPEEDSEPARATfOJJ 
HISTOGRAM  FROM  SINGLE  A fRCRAFT  DATA 

The  Conflict  Monitoring  Model  requires  the  probability  that  an 
aircraft  pair  has  a  crosstrack  closing  speed  and  a  crosstrack 
separation  within  a  given  range.  Since  the  FAA ' s  data 
collection  could  not  take  this  data  on  aircraft  pairs,  single 

aircraft  data  on  cros.strack  deviations  and  crosstrack  speed  that 
was  taken  was  used.  This  appendix  will  outline  the  procedure  to 
convert  a  joint  histogram  of  crosstrack  deviation  and  crosstrack 
speed  to  a  joint  histogram  of  crosstrack  separation  and 
crosstrack  closing  speed. 

First  let  us  assume  that  we  have  the  bivariate  probability 
density  function  of  the  crosstrack  deviation  and  crosstrack 

speed  from  single  aircraft  observations.  We  will  call  this 
probability  density  function  for  aircraft  on  route  i 

f[(Yj,Vy{l.  We  could  then  compute  another  probability 

density  function,  g,  a  function  of  y  =  Y2-Y},  Vyj,  and 

vy2  ai> 

o° 

g(y.Vyl,Vy2)  =  /  fjiYj.Vy!)  f2(y-Y1,Vy2)dY1  (F-l) 

-'h 

But  we  really  have  histogram  data  rather  than  a  probability 
density  function  so  that  we  will  denote  f{(Y^,Vy^)  as  Che 
estimate  for  f[(Y|,Vy|)  based  on  the  histogram  data 
f[(Yi,Vyi).  The  quantities  F  and  f  are  related  by  the 

following  expression: 

C  Yi , Vy f)=  fiCYi.Vyj)  (F-2) 

AY  AVy~ 

where  AY  is  the  histogram  cell  size  for  crosstrack  deviations 
and  AVy  is  the  histogram  cell  size  for  the  crosstrack  speed. 
We  can  then  write  the  estimate  for  the  function  g  as 
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But 
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Therefore 


( F-4 ) 
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The  expression  in  (P-5)  is  the  histogram  convolution  of  the 
crosstrack  deviations.  The  result,  g,  is  a  histogram  in  terms 
of  the  spacing  between  the  aircraft  and  their  individual 
crosstrack  velocities.  By  letting  y=  Vy2-Vyl,  we  can 
change  the  variables  in  the  function  g.  Thus  y 


2fy»vyl .Vy2)  =  g  ( v , Vy  2— y  ,Vy2) 


(F-h) 


By  summing  over  Vy2  we  arrive  at 

g(y,y'>  ’  2fy,  vy2-y,  vy2)  (f-7) 

''■’vj  =*-00 

The  histogram  g(y,y)  is  the  desired  histogram.  The  histogram  is 
centered  on  a  crosstrack  closing  speed  of  zero  and  a  crosstrack 
separation  of  zero.  To  use  the  histogram  it  ha3  to  be  shifted 
so  that  its  center  corresponds  to  the  spacing  between  the  routes 
to  which  the  aircraft  have  been  assigned.  Once  the  histogram 
g(y,$)  has  been  computed,  the  probability  that  an  aircraft  pair 
has  a  given  range  of  crosstrack  closing  speed  and  crosstrack 
separation  can  be  estimated  by  adding  up  the  appropriate  cells 
of  the  histogram. 
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I NTEKVKNT  ION  RATI-  S  IMU1.ATION 

Cl  .  INTRODUCTION 

This  appendix  will  discuss  tin*  estimation  of  the  controller 
intervention  rate.  The  est  imate  of  this  system  performance 
measure  was  obtained  through  a  simulation.  The  input  data  to 
the  simulation  will  be  described  followed  bv  an  outline  of  the 
steps  the  simulation  performs  to  arrive  aL  the  intervention  rat- 
estimates.  The  analysis  of  the  output  f  rum  the  simulation  is 
dealt  with  at  the  end  of  tin-  appendix. 

C.2  .  THE  INPUT 

The  data  which  were  used  in  the  simulation  were  smoothed, 

aircraft  tracks  which  are  part  of  the  FAA 1 s  VOR  navigation  data 
collection.  About  100  aircraft  on  each  of  two  selected  routes 
in  the  Cleveland  ARTCC  were  randomly  chosen  over  the  time  period 
of  the  data  collection.  The  only  requirement  on  these  aircraft 
were  that  they  had  actually  f ' own  the  entire  route  segments  of 
interest.  The  other  data  us-  i  in  the  simulation  are  shown  in 
Figure  G- 1 .  These  data  : n- lud-  the  radar  noise  and  the  location 
of  the  radar,  the  route  conr  ’nates,  the  sector  boundaries,  a 
random  number  seed,  the  number  of  route  spacings  and 

replications  for  the  simulation  to  consider,  and  the  requested 
flow  rates  of  the  traffic  on  the  routes.  The  time  duration 
the  simulation  is  also  spot  i  t‘  i  •*<.!. 

G3 .  THE  SIMULATION  FLOW 

The  object  of  the  simulation  is  to  estimate  the  controller 
intervention  rate  at  several  route  spacings  for  several  flow¬ 

rates  along  the  routes.  Since  this  is  a  simulation 
"experiment,"  several  replications  should  be  made  so  that 
meaningful  statements  an  be  made  concerning  the  results. 

To  minimize  the  cost  of  these  numerous  runs,  a  definite 

simulation  run  strategy  wa  •;  employed.  The  strategy  depends  on 
the  use  of  the  same  trait  ic  for  each  different  route  spacing 
during  the  same  replication  of  the  simulation.  The  minimum 
route  separation  to  be  considered  (8  nm i 1  was  run  first.  Of  the 
traffic  run  on  the  routes  spaced  a  minimum  distance  apart,  only 
a  subset  of  the  traffic  would  form  the  set  of  potential 

conflicts  for  the  next  set  of  wider  spaced  routes.  In  i  ther 

words,  with  the  same  traffic,  aircraft  which  are  not  projected 
to  get  near  to  other  aircraft  at  the  narrower  route  spacing  will 
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FIGURE  G-1 

INTERVENTION  RATE  SIMULATION  FLOW 
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not  be  projected  to  get  near  to  other  aircraft  when  the  routes 
are  farther  apart.  Of  course  some  care  must  be  taken  in 
choosing  the  definition  of  what  it  means  to  be  projected  to  get 
near  to  another  aircraft  for  the  purposes  of  carrying  the 
aircraft  along  to  the  next  route  spacing.  The  reason  for  this 
care  is  that  an  aircraft  pair  may  barely  not  be  detected  in 
conflict  at  one  route  spacing  but  might  be  detected  in  conflict 
at  a  wider  route  spacing.  For  instance,  at  the  narrow  route 
spacing,  a  deviating  aircraft  might  be  projected  to  pass  in 
front  of  an  aircraft  on  the  adjacent  route  but  with  enough 
separation  not  to  cause  an  alarm.  However,  with  the  aircraft  in 
exactly  the  same  positions  except  on  routes  that  are  spaced 
farther  apart,  the  projected  minimum  separation  might  be  small 
enough  to  trigger  an  alarm.  To  guard  against  not  including  such 
aircraft  in  the  wider  route  spacing  sample  an  additional 
"window"  which  is  larger  than  the  conflict  alert  criteria  was 
employed. 

The  flow  of  the  simulation  itself  is  rather  straightforward  as 
shown  in  Figure  G-l.  First,  the  traffic  is  generated.  This  is 
done  by  randomly  selecting  a  track  from  the  track  data.  This 
data  consists  of  the  best  estimate  of  an  aircraft's  position 
every  12  seconds  during  the  time  the  aircraft  was  observed  in 
the  FAA's  data  collection.  Next  an  entry  time  is  chosen  for 
that  aircraft.  The  entry  time  is  chosen  such  that  the 
interaircraft  entry  times  axe  exponentially  distributed 
corresponding  to  the  desited  average  flow  rate  on  the  route. 
Since  exponentially  distributed  interaircraft  entry  times  could 
allow  two  aircraft  to  be  very  closely  separated  in  time  along 
the  same  route,  a  check  is  made  to  insure  that  the  radar 
separation  is  not  violated  between  aircraft  on  the  same  route. 
If  an  aircraft  pair  is  detected  to  violate  the  separation 
standard,  the  second  aircraft  is  delayed  at  the  route  entry 
point  to  insure  the  minimum  separation,  plus  a  buffer,  between 
the  two  aircraft. 

After  the  traffic  is  generated,  the  coordinates  of  the  aircraft 
positions  are  transformed  to  correspond  to  two  routes  that  are 
spaced  at  the  minimum  spacing.  At  this  point  the  simulation  has 
a  string  of  positions  and  times  that  is  organized  by  aircraft. 
To  sismlate  radar  returns,  the  aircraft  positions  must  be 
organized  by  time.  This  is  achieved  by  sorting  the  aircraft 
position  records  by  time  (radar  update). 

The  simulation  can  now  track  and  perform  conflict  alert  on  the 
aircraft  positions.  But  before  the  simulation  can  perform 
either  of  these  functions  on  a  particular  set  of  aircraft. 
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certain  conditions  must  prevail.  The  first  condition  is  that 
the  time  of  the  radar  update  must  be  between  the  starting  and 
ending  tiae  values  of  the  simulation.  These  values  are  set  when 
the  traffic  is  generated  to  account  for  the  "end  effects"  of  the 
simulation.  The  simulation  does  not  start  until  the  latter  of 
the  first  aircraft  on  each  route  exits  the  sector.  The 
simulation  ends  when  the  earlier  of  the  last  aircraft  on  each 
route  enters  the  sector.  After  the  starting  update,  the  tracker 
will  maintain  a  track  on  all  the  aircraft. 

The  positions  of  the  aircraft  from  the  "track  data"  have  been 
saoothed  as  part  of  the  FAA's  data  collection  program. 
Therefore,  radar  noise  (azimuth  and  tange)  is  added  to  the 
position  reports  prior  to  going  into  the  NAS  tracker  in  the 
simulation.  The  tracked  aircraft  are  then  subjected  to  the 
conflict  alert  function  if  two  further  conditions  are  met.  One 
condition  is  that  both  aircraft  have  to  be  within  the  sector  of 
interest.  The  other  condition  is  that  there  has  not  been  a 
previous  conflict  alert  declared  on  this  aircraft  pair. 

An  aircraft  pair  which  satisfies  the  above  conditions  are  then 
passed  to  the  conflict  alert  function  where  the  pair  is 
subjected  to  a  set  of  coarse  geographic  and  velocity  filters. 
Passing  these  coarse  filters  indicates  that  the  pair  is  near  and 
generally  closing  on  each  other.  The  aircraft  pairs  that  pass 
the  coarse  filtering  are  then  subjected  to  a  set  of  fine 
filters.  These  fine  filters  project  the  positions  of  the  pair 
ahead  in  tiae.  If,  within  a  given  look-ahead  time,  the  pair  is 
projected  to  be  separated  by  less  than  a  certain  distance,  then 
the  pair  is  called  a  conflict  pair.  If  the  pair  passes  the  fine 
filtering  in  two  out  of  the  last  three  passes  through  the 
filters  then  the  pair  is  in  potential  conflict.  In  the  NAS,  the 
controller  is  alerted  to  this  conflict  pair  by  the  NAS 
automation  blinking  the  aircraft  symbology  on  his  screen.  In 
the  siaulation,  the  potential  conflict  is  tallied  as  part  of  the 
simulation  statistics.  The  current  conflict  and  predicted 
conflict  parasieters  are  the  same  ones  used  in  NAS.  In  addition 
other  parameters  are  used  by  the  simulation  to  identify  those 
aircraft  pairs  which  should  be  included  in  the  next  wider  route 
spacing  run.  The  definitive  formulation  of  the  NAS  Conflict 
Alert  function  can  be  found  in  Reference  (13).  Appendix  H  gives 
a  detailed  description  of  the  horizontal  components  of  the  NAS 
Conflict  Alert. 

Referring  back  to  Figure  G-l,  we  can  see  that  after  the  conflict 
alert  has  been  performed,  we  have  a  list  of  nearly  conflicting 
pairs  of  aircraft  to  be  used  in  the  simulation  of  the  next  wider 
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route  spacing.  If  there  are  no  more  nearly  conflicting  pairs  or 
there  are  no  more  route  spacings  to  consider,  we  go  back  to  do 
another  replication  of  the  simulation.  If  there  are  nearly 
conflicting  pairs  and  more  route  spacings  then  the  simulation 
reassigns,  the  nearly  conflicting  pairs  to  the  routes.  The  same 
entry  times  that  they  had  originally  are  maintained.  The  newly 
assigned  aircraft's  positions  are  then  transformed  to  correspond 
to  the  new  route  spacing.  The  loop  is  thus  closed  and  the 
processing  continues  until  all  the  route  spacings  and 
replications  are  done. 

G4 .  THE  OUTPUT 

The  output  from  the  simulation  consists  of  several  statistics. 
First,  the  traffic  generator  lists  the  traffic  characteristics: 
the  particular  aircraft  chosen,  their  entry  times  onto  the 
routes,  the  delays  that  were  incurred,  and  the  realized  average 
traffic  flow  rates.  The  tracking  and  conflict  alert  part  of  the 
simulation  records  another  set  of  statistics.  These  include 
the  number  of  flying  hours  in  the  sector,  the  number  of  aircraft 
generating  those  flying  hours,  the  maximum  instantaneous 
aircraft  count  (IAC),  the  number  of  hours  the  conflict  alert 
function  was  "watching"  the  sector,  and  the  number  of  conflict 
alerts  declared.  Other  miscellaneous  statistics  concerning  the 
detailed  workings  of  conflict  alert  are  also  taken. 

G5 .  THE  ANAYSIS  OF  THE  OUTPUT 

As  Reference  (14)  points  out,  there  are  basically  two  types  of 
simulations  as  far  as  the  analysis  of  the  simulation  output  is 
concerned  —  terminating  and  steady-state.  The  difference 
between  these  is  the  desired  result.  In  all  cases  a  simulation 
has  to  be  terminated  after  some  finite  time  interval.  In  some 
cases,  one  wishes  to  evaluate  steady  state  property  of  the 
system  as  if  the  time  interval  goes  to  an  infinite  length.  In 
other  cases,  the  behavior  of  the  system  over  a  defined 
(terminated)  length  of  time  is  desired  such  as  the  peak  of  the 
traffic  demand.  The  simulation  of  conflict  alerts  belongs  to 
the  terminating  type  of  simulation.  One  of  the  results  of 
Reference  (14)  is  that  for  a  terminating  simulation  the  usual 
statistical  tests  and  procedures  can  be  applied  to  the  output 
from  the  simulation. 

In  the  particular  simulation  described  above,  the  event  of  a 
conflict  alert  occurs  over  a  period  of  time.  It  is  also  as 
likely  that  this  event  will  occur  in  one  interval  as  in  any 
other  and  that  the  occurrence  of  the  event  has  no  effect  on 
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whether  or  not  another  occurs.  The  number  of  events  in  a  fixed 
time  with  the  above  attributes  is  often  assumed  to  have  a 
Poisson  distribution.  In  fact  in  the  results  section  (Section  5 
in  Volume  I  of  this  report)  it  is  shown  that  the  number  of 
conflict  alerts  per  hour  from  the  simulation  does  look  as  if  it 
were  Poisson  distributed. 

The  result  we  are  looking  for  is  the  expected  value  of  the 
Poisson  distribution.  The  unbiased  and  maximua  likelihood 
estimator  is  the  sample  mean.  This  is  the  number  of  conflict 
alerts  divided  by  the  number  of  hours  which  were  simulated.  We 
can  also  construct  a  confidence  interval  for  the  expected  value 
of  a  Poisson  distribution  (Reference  (15).  If 
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f(n)  *  — j  n  -  0,1,2... 


(G-l) 


then  for  any  value  n1  and  o<.5,  lower  and  upper  limits  of  may 
be  determined  such  that 
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such  that  <  Xj,.  The  values  of  and  are  tabulated 
for  l00(l-2»)%  »  95%  and  99%  significance  levels  in  Reference 
(15).  These  tables  are  applied  to  the  simulation  output  to 
construct  95%  confidence  intervals  for  the  expected  rate  of 
conflict  alerts.  The  results  are  shown  in  Section  5  of  Volume  I 
of  this  report. 
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APPENDIX  H 


THE  NAS  CONFLICT  ALERT 


In  Section  3.1.1  of  Volume  I  of  this  report  it  was  stated  that 
the  conflict  alert  function  in  the  NAS  computer  was  not  modelled 
precisely  in  the  Conflict  Monitoring  Analysis.  The  NAS  Conflict 
Alert  function  uses  both  horizontal  and  vertical  filters  when 
determining  the  condition  of  conflict  between  a  pair  of 
aircraft.  Since  we  have  assumed  the  aircraft  to  be  in  straight 
and  level  flight  at  the  same  altitude,  this  appendix  will 
discuss  only  the  horizontal  filters  used  by  the  NAS  Conflict 
Alert  function. 

In  the  NAS  computer  the  Conflict  Alert  performs  a  series  of 
tests  to  determine  the  conflict  status  of  a  pair  of  aircraft. 
The  tests  are  structured  in  such  a  way  as  to  minimize  false 
alarms,  maximize  the  chance  of  an  alarm  on  a  true  conflict  while 
at  the  same  time  minimizing  the  computational  complexity  of  the 
tests. 

After  the  tracks  are  geographically  sorted  and  tested  for  being 
near  in  altitude,  a  series  of  horizontal  filters  are  applied. 
The  filters  are  diagramed  in  Figure  H-l.  The  Conflict  Alert 
has  two  types  of  horizontal  filters,  coarse  and  fine.  For  each 
type  there  may  be  two  levels  designated  as  A  and  B.  The  type 
and  level  of  the  filter  will  depend  on  the  previous  history  of 
the  aircraft  pair  with  respect  to  the  Conflict  Alert  function. 
The  values  of  the  Conflict  Alert  parameters  are  summarized  on 
Table  H-l. 

The  process  starts  with  position  estimates  (4X£,  Xj,  Y^ , 
Y;)  and  velocity  component  estimates  Xj,  Yj,  Yj) 

of  aircraft  i  and  j.  It  is  first  determined  what  type  and  level 
of  filter  should  be  used.  The  choice  will  depend  on  whether  the 
aircraft  pair  has  been  observed  previously  and,  if  so,  what  the 
result  of  the  previous  Conflict  Alert  tests  were.  If  the 
aircraft  pair  has  not  passed  any  filters  successfully  in  the 
past  two  applications  of  the  filters,  a  level  A  filter  is  used. 
If  a  level  A  filter  is  required  then  a  coarse  horizontal  filter 
is  performed  on  the  position  estimates: 

(Xi-Xj)2  +  (Yi-Yj)2<  MAXR2  (H-l) 

where  MAXR“55  nmi.  If  the  aircraft  are  separated  by  more  than 
5$  nmi  they  are  no  longer  considered  by  Conflict  Alert. 


FIGURE  H-1 

HORIZONTAL  CONFLICT  ALERT  FILTERS  i 
USED  IN  SIMULATION  I 


HORIZONTAL  CONFLICT  ALERT  PAPA BETEFS 


If  Che  coarse  filter  is  passed  the  positions  and  velocities  of 
the  aircraft  pair  are  subjected  to  the  fine  filter  of  level  A. 


The  fine  filter  first  checks  if  the  aircraft  pair  is  in  current 
lateral  conflict.  If 

Rg  -  (XrXj)2  ♦  (Yi-Yj)2  (H-2) 

is  less  than  SEPR2,  a  current  horizontal  conflict  exists.  For 
level  A  SEPR  *  5  nmi  and  for  level  B  SEPR  =4.2  nmi. 

If  the  aircraft  pair  is  not  in  current  conflict,  it  could  be  in 
predicted  conflict.  This  part  of  the  fine  filter  projects  an 
aircraft  pair  ahead  in  a  straight  line  for  a  certain  period  in 
time  and  tests  to  see  if  the  aircraft  come  within  a  given 
distance  of  each  other  during  that  time  interval.  This  check  is 
done  in  several  steps  to  reduce  computation  time.  The  first 
step  is  to  determine  whether  the  tracks  are  generally  converging 
toward  each  other.  If  we  let 

V^CXi-XjHii-ijMYi-YjXVtj)  (H-3) 

then  the  inquality  Vc  <  VELC  must  be  satisfied  to  continue 

testing.  For  both  level  A  and  level  B,  VEL08  nmi2/hr.  The 
sign  of  quantity  Vc  indicates  a  closing  condition  in  the  X 
and  Y  directions.  If  the  aircraft  were  closing  in  both  X  and  Y 
simultaneously  then  tc  would  be  less  than  zero.  If  the 
aircraft  were  opening  in  both  X  and  Y  simultaneously,  then  Vc 
would  be  positive.  If  the  aircraft  were  closing  in  one 

direction  and  opening  in  the  other  direction,  the  value  of  Vc 
may  be  positive  or  negative.  Even  if  Vc  is  positive  the  miss 
distance  might  be  small  enough  to  qualify  as  a  potential 
conflict.  Thus  Vc  is  tested  against  a  positive  nimtber. 

If  the  aircraft  are  closing,  they  may  be  doing  so  very  slowly. 
If  the  closure  is  too  slow  there  cannot  be  a  conflict  within  the 
prescribed  time  interval.  If  we  let 

v2*(xrXj)2  +  (frtj)2  (H-4) 

then  the  inquality  v2  >  CL0S2  must  be  satisfied  to  continue 

testing.  CLOS  is  equal  to  .7  nmi/min  for  both  level  A  and  level 
B. 

The  predicted  minimum  separation  is  tested  next.  The  time  to 
minimum  separation,  Tm,  is  -Vc/v2.  if  v2*0,  then 
Tra*0.  The  square  of  the  predicted  minimum  separation  is 
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If  t SEPM2  then  the  checking  can  continue.  For  level  A 
SEW1  is  set  to  6  nmi.  For  level  B  SEPM  is  set  to  4.8  nmi .  The 
window  is  opened  wider  on  level  A  to  allow  for  a  turning 
aircraft  to  be  detected  at  the  earliest  moment. 

After  it  is  determined  that  the  aircraft  are  predicted  to  close 
to  less  than  the  preset  minimum,  a  test  is  performed  to 
determine  if  the  minimum  separation  is  violated  within  the 
look-ahead  time  (WRNT).  For  level  A  WRNT  *  2.5  min.  For  level 
B,  WRNT  *  2.0  min.  The  minimum  separation  is  violated  within 
the  look-ahead  time  if  one  of  the  following  conditions  holds: 

•  The  time  to  minimum  separation,  Tm,  is  less  “than  or 
equal  to  WRNT  or 

•  The  square  of  the  separation  predicted  at  WRNT  violates 
the  separation  limits. 

The  square  of  the  separation  predicted  at  WRNT  is  given  by 
R£(WRNT)=R§+(2VC*WRNT)+V2*WRNT2  (H-6) 

Violation  of  the  separation  limits  occurs  when 
R2(WRNT)<  SEPP2. 

The  value  of  SEPP  for  level  A  is  6  nmi.  For  level  B  the  value 
is  4.8  nmi.  Here  again,  the  longer  time  window  on  level  A  is 
designed  to  allow  more  efficient  detection  of  accelerating 
aircraft. 

If  an  aircraft  pair  passes  both  the  fine  horizontal  filter  and 
the  fine  altitude  filter  (which  is  not  being  considered  here)  at 
level  A  then  it  is  placed  in  a  conflict  pairs  table.  On  the 
next  tracking  cycle  this  conflict  pair  will  be  tested  against 
the  level  B  set  of  filters.  An  aircraft  pair  is  eligible  for 
controller  alert  generation  if  it  has  been  determined  to  be  in  a 
condition  of  conflict  at  least  twice  in  the  past  three  sucessive 
applications  of  the  filters.  This  procedure  reduces  the  rate  of 
false  alarms  due  to  noise  in  the  surveillance  system.  Once  an 
aircraft  pair  has  been  in  a  condition  of  conflict  for  the 
requisite  number  of  successful  tests,  the  pair  will  be 
identified  to  the  controller  by  flashing  data  blocks  and  list 
displays  at  the  controller's  position.  This  flashing  will 
continue  as  long  as  the  aircraft  pair  continues  to  meet  the 
Conflict  Alert  criteria  or  the  controller  manually  turns  off  the 
alarm. 


The  parameters  in  Table  H~  1  labeled  as  Near  Conflict  parameters 
are  those  used  by  the  intervention  rate  simulation  described  in 
Appendix  G.  These  parameters  are  used  to  identify  those  pairs 
°f  **rC.r*^C  ^at  might  he  potential  conflict  pairs  in  the 
simulation  run  at  a  wider  route  spacing. 
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The  tot.il  delay  time  at  which  the  aircraft 
pair  will  come  into  horizontal  overlap 
before  there  is  an  avoidance  turn.  (Figure 
3-5 ,  Vo  1  ume  I ) . 

One  corner  of  the  polygon  B  (Figure  C-4, 
Volume  IT). 

The  alongtrack  positions  of  two  aircraft  on 
their  respective  routes  (Figure  E-l,  Volume 
IT). 

Azimuth  Count  Pulse.  A  unit  of  angular 
measure  equal  to  1/4096  of  a  circle  (.0379 
degrees)  (Table  5—1,  Volume  T). 

A  parameter  of  the  integral  of  the  circular 
normal  probability  function  (Equation  C-6, 
Vol ume  IT). 

The  polygonal  area  over  which  one  computes 
the  probability  of  observing  an  aircraft 
pair  (Figure  C-4 ,  Volume  TT). 

One  corner  of  the  polygon  B  (Figure  C-4, 

Volume  TT). 

The  alongtrack  position  of  an  aircraft  on  a 
route  (Figure  E-l,  Volume  IT). 

The  starting  time  of  the  periodic 
observations  (Section  D-4,  Volume  II). 

A  parameter  used  to  fit  the  envelope  for  the 
horizontal  overlap  region  (Equation  B-21, 
Volume  TT) 

Closing  Speed  Threshold  parameter  (Equation 
H-4 ,  Volume  IT). 

One  corner  of  the  polygon  B  (Figure  C-4, 

Volume  IT). 


The  minimum  projected  separation  between  an 
aircraft  pair  u9ed  by  the  conflict  alert 
function  (Equation  3-4,  Volume  I). 

One  corner  of  the  polygon  B  (Figure  C-4 , 
Volume  TI). 

The  separation  between  two  aircraft 
(Equation  A-l,  Volume  II). 

The  minimum  separation  between  a  pair  of 
aircraft  (Figure  3-3,  Volume  I). 

The  expected  value  of  the  separation  between 
aircraft  on  the  same  route  (Equation  E-3, 
Volume  IT). 

The  crosstrack  distance  inside  the  conflict 
region  at  the  initial  radar  observation 
(Figure  B-2,  Volume  II).  The  partitioning 
parameter  for  the  crosstrack  closing  speed 
(Equation  B-18,  Volume  II). 

The  Fast  Fourier  Transform  (Appendix  D, 
Volume  IT). 

The  Poisson  probability  density  function 
(Equation  G-l,  Volume  II). 

A  finite  Fourier  transform  (Equation  (D-l) , 
Volume  IT. 

A  finite,  complex  series  (Equation  D-l, 
Volume  II). 

The  probability  density  function  of  the 
separation  between  two  aircraft  on  the  same 
route  (Equation  E-l,  Volume  II). 

The  histogram  data  corresponding  to  “fj . 
(Equation  F-2,  Volume  II). 

The  even  elements  of  the  transformed  series, 
F(n)  (Equation  D-6,  Volume  II). 

The  even  elements  of  the  series  f(t) 
(Equation  D-6,  Volume  II). 
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The  joint  probability  density  function  of 
crosstrack  deviation  and  crosstrack  speed 
from  single  aircraft  observations  from 
aircraft  on  route  i  (Equation  F-l,  Volume 
TT). 

The  estimate  for  f;(Y;  Vvj)  (Equation 
F-2,  Volume  IT).  ’  y 

The  histogram  data  associated  with  the 
probability  density  function 

fi(Vi,Vyi)  (Equation  F-2,  Volume  IT). 

The  odd  elements  of  the  transformed  series, 
F(n)  (Equation  D-6,  Volume  II). 

The  odd  elements  of  the  series  f(t) 
(Equation  D-b,  Volume  II). 

The  crosstrack  separation  lost  during  one 
radar  update  time  interval  (Figure  B-2, 
Volume  II). 

The  gravitational  constant  (Equation  B-6, 
Volume  II). 

A  finite  Fourier  Transform  (Equation  D-5, 
Volume  II). 

A  finite,  complex  series  (Equation  D-5, 
Volume  IT). 


g(y,y) 


S(y»vyl*Vy2) 


sfy>vyi»vy2) 


The  derived  histogram  of  crosstrack 
separation  and  crosstrack  closing  speed. 
(Equation  F-7,  Volume  IT). 

The  joint  probability  density  function  of 
the  cross  track  separation  and  the  two 
crosstrack  speeds  (Equation  F-l,  Volume  II). 

The  estimate  for  g(y,Vvl,Vv2)  (Equation 
F-3 ,  Volume  II).  r  * 


?(y,Vyi,Vy2) 


The  histogram  data  corresponding  to  £. 
(Equation  F-4,  Volume  II). 


Dummy  Variable  (Equation  B-14,  Volume  II). 
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A  distance  defined  in  Figure  C-6,  Volume  It 

The  convolved  density  function  (Equation 
D-ll,  Volume  II). 

The  probability  density  function  of  the  bank 
angle  k  (Figure  B-4,  Volume  II). 

The  probability  density  function  of  the 
alongtrack  distance  between  aircraft  on  one 
route  and  an  aircraft  on  an  adjacent  route 
(Equation  E-2,  Volume  II). 

Dummy  Variable  (Equation  B-14,  Volume  II). 

Instantaneous  aircraft  count  (Appendix  G, 
Volume  II). 

Dummy  Variable  (Equation  B-14,  Volume  II). 

The  probability  density  function  of  the  turn 
rate  ^(Equation  B-7,  Volume  II). 

Dummy  Variable  (Equation  B-14,  Volume  It). 

Distances  defined  in  Figure  C-6,  Volume  IT. 

i 

Dummy  Variable  (Equation  B-14,  Volume  II). 

Dummy  Variable  (Equation  B-14,  Volume  II), 
(Appendix  A,  Volume  II). 

The  maximum  separation  at  which  Conflict 
Alert  will  continue  to  consider  a  pair  of 
aircraft  (Equation  H-l,  Volume  II). 

The  half  series  length  (Equation  D-6,  Volume 
II),  also  the  number  of  equally  spaced 
initial  positions  in  the  Conflict  Region 
Boundary  (Equation  B-5,  Volume  II). 

Index  of  the  FFT  (Equation  D-l,  Volume  II). 

The  average  flow  rate  of  aircraft  on  routes 
1  and  2  (Equation  3-7,  Volume  I). 

The  factor  which  delineates  the  three  cases 
of  overlap  (Equation  B-ll,  Volume  II). 
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The  probability  of  first  detecting  an 

aircraft  pair  on  update  i.  (Equation  B-4 , 
Volume  IT). 

The  probability  of  first  detecting  an 

aircraft  pair  on  update  i  given  the  first 
observation  was  at  F.  (Equation  B-2,  Volume 
IT). 

The  proportion  of  time  that  the  horizontal 
separation  is  lost  (Section  2.2,  Volume  I). 

The  probability  that  the  time  delay  has  a 

value  in  the  kth  cell.  This  probability 
will  also  depend  on  the  conflict  region 
boundary  cell,  i.  (Equation  B-27,  Volume  II). 

The  probability  that  the  turn  rate  has  a 

value  in  the  £th  cell  (Equation  B-27,  Volume 
II). 

The  probability  that  the  aircraft  pair  will 
not  be  detected  in  the  conflict  region  on 
radar  update  i  given  the  first  observation 
was  at  F  (Equation  B.-J  ,  Volume  II). 

The  probability  of  not  being  detected  in  the 
conflict  region  during  the  first  i  updates 
given  the  first  observation  was  at  F. 
(Equation  B-l,  Volume  II). 

The  proportion  of  time  that  the  alongtrack 

separation  is  lost  (Section  2.1,  Volume  I). 

The  proportion  of  time  that  the  crosstrack 

separation  is  lost  (Section  2.1,  Volume  I). 

The  probability  of  an  aircraft  pair  being  on 
the  conflict  region  boundary  (Equation  B-23, 
Volume  IT). 

The  probability  of  an  aircraft  pair  being  on 
the  conflict  region  boundary  in  integration 
cell  i.  (Equation  3-2,  Volume  I). 

The  probability  of  an  aircraft  pair  coming 
into  horizontal  overlap  given  that  the  pair 
started  on  the  conflict  region  boundary  in 
integration  cell  i  (Equation  3-2,  Volume  I). 
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p<Yi  »Ayi *Yi  *A$^ 

The  joinf  probability  that  the  aircraft  pair 
has  a  crosstrack  separation  y^  in  the  i  th 
cell  and  has  a  crosstrack  closing  speed  yj 
in  the  ith  celt  (Equation  B-26,  Volume  III. 

PA< 

The  probability  that  an  aircraft  pair  has  an 
interroute  alongtrack  separation  Ax 
(Equation  3-4,  Volume  I). 

% 

The  number  of  multiplications  required  for 
calculating  the  transform  of  air  N  point 
array.  ^Equation  D-7,  Volume  II). 

q(F) 

The  probability  density  function  of  the 
distance  inside  the  conflict  region  boundary 
at  which  the  first  observation  is  made 
(Equation  B-3,  Volume  IT). 

R 

The  radius  of  the  right  cylindrical 

collision  shape  (Section  3.2,  Volume  T). 

Rm 

The  predicted  minimum  Separation  (Equation 
H-S,  Volume  II). 

Ro 

The  lateral  separation  between  tracks 

(Equation  H-2,  Volume  II). 

RP 

Predicted  separation  at  the  warning  time 
(Equation  H-6,  Volume  II). 

rv2i 

The  proportion  of  the  k  th  cell  in  the 
overlap  space  is  enclosed  within  the 

horizontal  overlap  region  corresponding  to 
the  ith  conflict  region  boundary  cell 
(Equation  B-27,  Volume  II). 

SEPM 

Minimum  Separation  Parameter  (Equation  H-5, 
Volume  II). 

SEPP 

The  separation  at  the  warning  time  parameter 
(Equation  H-6,  Volume  II). 

SEPR 

The  current  lateral  conflict  threshold 

(Equation  H-2,  Volume  II). 

T 

The  number  of  elements  in  the  series  on 
which  the  FFT  operates  (Equation  D-l,  Volume 
II). 

The  time  variable,  also  a  summation  index 

The  total  delay  time  from  the  time  the 
aircraft  pair  enters  the  conflict  region  to 
the  time  one  aircraft  starts  its  avoidance 
maneuver.  (Equation  B-10,  Volume  IT). 

The  look-ahead  time  in  the  conflict  alert 
function  (Equation  A-3,  Volume  IT). 

Time  to  minimum  separation  (Equation  H-5, 
Volume  II). 

The  time  the  aircraft  ia  in  the  turn  to  the 
point  of  closest  approach  (Equation  B-10, 
Volume  IT). 

The  time  the  aircraft  pair  is  at  the 
conflict  region  boundary  (Figure  3-3,  Volume 
I). 

The  time  at  which  the  system  has  cognizance 
that  the  aircraft  pair  is  within  the 
conflict  region  (Figure  3-3,  Volume  I). 

The  time  at  which  the  aircraft  start  its 
avoidance  maneuver  (Figure  3-3,  Volume  I). 

A  coordinate  axis  used  to  circularize  the 
bivariate  normal  distribution  (Equation  C-3, 
Volume  II). 


The  u  coordinate  value  of  point  A'  (Equation 
C-9,  Volume  II). 

The  u  coordinate  value  of  point  B’  (Equation 
C-9,  Volume  II). 

The  forward  velocity  of  the  aircraft 
(Equation  3-7,  Volume  I). 

A  coordinate  axis  used  to  circularize  the 
bivariate  normal  distribution  (Equation  C-3, 
Volume  IT). 

Converging  speed  indicator  (Equation  H-3, 
Volume  II). 

Converging  speed  threshold  parameter 
(Equation  H-3,  Volume  II). 
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Closing  speed  indicator  (Equation  H-4, 
Volume  II). 

The  tracker  x  velocity  estimate  for  scan  n 
(Equation  C-l,  Volume  IT). 

The  alongtrack  speeds  of  aircraft  on  routes 
1  and  2  (Figure  A-l,  Volume  II). 

The  maximum  croastrack  speed  observed  in  the 
data  (Equation  A-4,  Volume  II). 

The  tracker  Y  velocity  estimate  for  scan  n 
(Equation  C-l,  Volume  II). 

The  crosstrack  speeds  of  aircraft  on  routes 
1  and  2  (Figure  A-l,  Volume  II). 

t 

The  integral  of  the  circular  normal 
probability  function  (Equation  C-6,  Volume 
II). 

The  v  coordinate  value  of  point  A'  (Equation 
C-9,  Volume  IT). 

The  v  coordinate  value  of  point  B'  (Equation 
C-9,  Volume  11). 

The  separation  between  two  aircraft  on  the 
same  route  (Figure  E-l,  Volume  II). 

Warning  Time  (Look-ahead  time)  (Equation 
H-A,  Volume  II). 

The  alongtrack,  interroute,  separation 
between  a  pair  of  aircraft  (Section  3.2, 
Volume  I). 

The  alongtrack,  interroute,  closing  speed 
between  a  pair  of  aircraft  (Equation  A-l, 
Volume  It). 

The  bounds  on  the  value  x  (Equation  A-7, 
Volume  II). 

The  alongtrack  position  of  an  aircraft  on 
route  i  (Figure  A-l,  Volume  II). 
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The  alongtrack  component  speed  of  an 
aircraft  on  route  i  (Equation  H-3,  Volume 
II). 

The  alongtrack  tragectory  of  the  aircraft  on 
route  i  (Equation  B-9,  Volume  II). 

The  alongtrack,  interroute,  separation  of  an 
aircraft  pair  in  the  ith  integration  cell 
(Figure  3-4,  Volume  I). 

The  tracker  predicted  X  position  estimate 
made  at  scan  n  for  the  position  at  scan  n*l 
(Equation  C-l,  Volume  II). 

The  reported  X  position  to  the  tracker  from 
the  surveillance  sjrstem  t  scan  n  (Equation 
C-l,  Volume  II). 

Ihe  tracker  X  position  estimate  for  scan  n 
(Equation  C-l,  Volume  II). 

The  X  and  Y  positions  of  an  aircraft  on 
route  i  when  it  is  on  the  conflict  region 
boundary  in  Conjunction  with  another 
aircraft  on  the  other  route.  (Equation  B-9, 
Volume  II). 

A  series  of  data  points  (Equation  D-8, 
Volume  II). 

The  crosstrack  separation  between  a  pair  of 
aircraft  (Section  3.2,  Volume  I). 

The  crosstrack  closing  speed  between  a  pair 
of  aircraft  (Section  3.2,  Volume  I). 

The  crosstrack  deviation  of  an  aircraft  from 
the  centerline  of  route  i  (Figure  A-l, 
Volume  II).  Also  the  Y  component 
(crosstrack)  position  of  an  aircraft  on 
route  i  (Equation  H-l,  Volume  II). 

The  crosstrack  component  speed  of  an 
aircraft  on  route  i  (Equation  H-3,  Volume 
II). 
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Yi(t) 


The  crosstrack  trajectory  of  the  aircraft  on 
route  i  (Equation  B-9,  Volume  II). 
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The  crosstrack  separation  of  an  aircraft 
pair  in  the  ith  integration  cell  (Figure 

3-4,  Volume  I). 

The  crosstrack  closing  speed  of  an  aircraft 
pair  in  the  ith  integration  cell  (Figure 

3-4,  Volume  I). 

The  crosstrack  separation  at  radar  update  i 
after  crossing  into  conflict  region  and 
being  initially  observed  at  F  (Figure  B-2, 
Volume  II). 

The  tracker  predicted  Y  position  estimate 
made  t  acan  n  for  the  position  at  scan  n+1 
(Equation  C-l,  Volume  II). 

The  reported  Y  position  to  the  tracker  from 
the  surveillance  system  at  scan  n  (Equation 
C-l,  Volume  II). 

The  tracker  Y  position  estimate  from  scan  n 
(Equation  C-l,  Volume  II). 

The  crosstrack  separation  at  the  conflict 
region  boundary  (Figure  B-2,  Volume  II). 

A  shifted  series  of  data  points  (Equation 
D-8,  Volume  II). 

A  variable  which  is  minimised  and  maximized 
in  Appendix  A.  (Equation  A-6,  Volume  II). 

The  convolved  density  function  (Equation 
D-ll ,  Volume  II). 

The  NAS  tracker  position  gain  (Equation  C-l, 
Volume  II),  and  the  significance  level  for 
he  confidence  interval  of  the  expected  value 
of  he  Foisson  distribution  (Equation  G-2, 
Volume  II). 

The  NAS  tracker  velocity  gain  (Equation  C-l, 
Volume  II). 

A  fixed  sampling  interval  (Section  D-4, 
Volume  II). 
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The  histogram  cell  size  for  crosstrack  speed 
(Equation  F-2,  Volume  II). 

The  integration  interval  in  the  nlongtrark 

dimension  (Figure  3-4,  Volume  I). 

The  integration  inverval  in  the  crosstrack 

dimension  (Figure  3-4,  Volume  I). 

The  integration  interval  in  the  crosstrack 

closing  speed  dimension  (Figure  3-4,  Volume 

1). 

The  histogram  cell  size  for  crosstrack 
deviations  (Equation  F-2,  Volume  IT). 

The  hank  angle  of  the  aircraft  making  the 
avoidance  maneuver  (Equation  B-6,  Volume  II). 

The  lower,  KL,  and  upper,  KU,  hounds  on  the 
hank  angle  distribution  (Equation  B-7, 

Vo 1 ume  T  T ) . 

A  parameter  of  the  integral  of  the  circular 
normal  probability  function  (Equation  C-6, 
Volume  T  T ) ,  and  the  parameter  of  the  Poisson 
distribution  (Equation  G-l,  Volume  II). 

The  lower  and  upper  confidence  limits  on  the 
expected  value  of  a  Poisson  distribution 
(Equation  G-2,  Volume  IT). 

The  length  of  the  rectangular  collision 
shape  which  represents  the  aircraft  (Table 
5-1,  Volume  T). 

The  width  of  the  rectangular  collision  shape 
which  represents  the  aircraft  (Table  5-1,  Volume 
I). 


The  expected  value  of  the  crosstrack  separation 
of  an  aircraft  pair  (Equation  C-2,  Volume  II). 

The  expected  value  of  the  crosstrack  closing 
speed  of  an  aircraft  pair  (Equation  C-2,  Volume 
II). 
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The  .ingle  the  velocity  vector  of  the  aircraft  on 
route  2  makes  with  the  crosstrack  axis  (Figure 
B-6,  Volume  II). 

The  angle  the  velocity  vector  of  the  aircraft  on 
route  1  makes  with  the  crosstrack  axis  (Figure 
B-6,  Volume  II). 

The  correlation  coefficient  between  the  errors 
in  the  crosstrack  closing  speed  and  separation 
(Equation  C-2,  Volume  II). 

The  normal  probability  integral  (Equation  C-8, 
Volume  II). 

The  turn  rate  of  the  aircraft  making  the 
avqidance  turn  (Equation  B-6,  Volume  II). 


APPENDIX  J 


REFERENCES 

1.  Flener,  Wiliam  M.  ,  "Request  for  E&D  Effort,"  Letter  from  the  l'AA 
Associate  Administrator  for  Air  Traffic  and  Airway  Facilities  to  the  FAA 
Associate  Administrator  for  Engineering  and  Development,  November  9, 
1976. 

2.  Kirkendall,  N.J.,  "A  Specification  for  a  Data  Collection  to 

Determine  Lateral  Pathkeeping  of  Aircraft  Flying  VOR  Defined  Routes  in 
High  Altitude  En  Route  Areas,"  The  MITRE  Corporation,  MTR-7430, 

Washington,  D.C.,  June  1977. 

3.  Kirkendall,  N.J.,  and  Smith,  A.P.,  "An  Extension  to  the  Data 

Specification  for  the  High  Altitude  Lateral  Pathkeeping  Study,"  The 
MITRE  Corporation,  MTR-7638,  Washington,  D.C.,  December  1977. 

4.  Colamosca,  B.F.,  Arch,  H.L. ,  Parkins,  K.D.,  "Data  Collection  for 
Study  of  Lateral  Separation  of  Aircraft  Flying  CONUS  High  Altitude 
VOR-Defined  Air  Routes,"  NAFEC  Technical  Letter  Report  NA-78-24-LR, 
RTCA  Paper  No.  69-78/SSRG-18 ,  May  1978. 

5.  Busch,  A.C.  and  Colamosca,  B.F.,  "Modeling  Safety  in  a  Procedural 

Air  Traffic  Control  Environment,"  Draft  NAFEC  Letter  Report 

NA-77-57-LR,  RTCA  Paper  No,  155-77/SSRG-l 1 ,  October  1977. 

6.  Polhemus,  N.W. ,  "Introduction  to  a  General  Collision  Risk  Model  for 
Intersecting  and  Nonparallel  Routes,"  RTCA  Paper  No.  9-79/SSRG-27 , 
January  1979. 

7.  George,  P.H.,  Johnson,  A.E.,  and  Hopkin,  V.D.,  "Radar  Monitorng  of 
Parallel  Tracks  —  Automatic  Warning  to  Controllers  of  Track  Deviations 
in  a  Parallel  Track  System,"  Eurocontrol  Experimental  Center  Report 
No. 67  (Task  C21/1),  September  1973. 

8.  Johnson,  N.L.,  and  Kotz,  S.,  "Distributions  in  Statistics: 
Continuous  Multivariate  Distributions,"  New  York:  John  Wiley  &  Sons, 
Inc.,  1972. 

9.  IEEE  Transactions  on  Audio  and  Electroacoustics,  Vol.  AU-15,  No.  2., 
June  1967. 

10.  Kleinrock,  L.,  "Queueing  Systems,  Volume  1:  Theory,"  New  York:  John 
Wiley  &  Sons,  Inc.,  1976. 


J-l 


11.  "Methodology  for  the  Derivation  of  Separation  Minima  Applied  to 
the  Spacing  Between  Parallel  Tracks  in  ATS  Route  Structures,  ICAO 
Circular  120-AN/89/2,  1976. 

12.  Kirkendall,  N.J.,  "Probability  of  Passings  Between  Aircraft  on 
Adjacent  Parallel  Routes:  Observed  Data  vs  Model  Estimates,"  The 
MITRE  Corporation,  WP-79W00397  Washington,  D.C.,  June  1979. 

13.  Hauser,  S.J.,  Dodge,  P.0,  and  Steinbacher,  J.G. ,  "Computer 
Program  Functional  Specifications  for  Conflict  Alert,"  The  MITRE 
Corporation,  WTR-7311,  Washington,  D.C.,  September  1976. 

14.  Law,  A.M. ,  "Statistical  Analysis  of  the  Output  Data  from 
Terminating  Simulations,"  Unpublished  and  Undated  Report  supported 
by  the  Office  of  Naval  Research  and  the  Army  Research  Office. 

15.  Handbook  of  Tables  for  Probability  and  Statistics.  The  Chemical 
Rubber  Co.,  Cleveland,  Ohio,  1968. 


J-2 


